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Space-time paraproducts for paracontrolled 
calculus, 3d-PAM and multiplicative Burgers 

equations 


I. BAILLEUlE, F. BERNICOlHand D. FReMI 


Abstract. We sharpen in this work the tools of paracontrolled calculus in order to 
provide a complete analysis of the parabolic Anderson model equation and Burgers 
system with multiplicative noise, in a 3-dimensional Riemannian setting, in either 
bounded or unbounded domains. With that aim in mind, we introduce a pair of 
intertwined space-time paraproducts on parabolic Holder spaces, with good continuity, 
that happens to be pivotal and provides one of the building blocks of higher order 
paracontrolled calculus. 
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Introduction 

It is probably understated to say that the work [26] of Hairer has opened a new 
era in the study of stochastic singular parabolic partial differential equations. It 
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provides a setting where one can make sense of a product of a distribution with 
parabolic non-positive Holder regularity index, say a, with a function with non¬ 
negative regularity index, say b, even in the case where o -I- 6 is non-positive, and 
where one can make sense of, and solve, a large class of parabolic stochastic singular 
partial differential equations by fixed point methods. The parabolic Anderson model 
equation (PAM) 

{dt + L)u = uC, (1.1) 

studied in Section [5] in a 3-dimensional unbounded background, is an example of 
such an equation, as it makes sense in that setting to work with a distribution 
of Holder exponent a — 2, for some a < ^, while one expects the solution u to the 
equation to be of parabolic Holder regularity a, making the product uC ill-defined 
since a -I- (a — 2) ^ 0. 

The way out of this quandary found by Hairer has its roots in Lyons’ theory of 
rough paths, which already faced the same problem. Lyons’ theory addresses the 
question of making sense of, and solving, controlled differential equations 

dzt = Vi{zt) dXi (1.2) 

in say, driven by an M^-valued ^-Holder control X = [X^,... ,X^'), with p ^ 2, 

and where V) are sufficiently regular vector fields on M'’*. Typical realizations of a 
Brownian path are ^-Holder continuous, with p > 2, for instance. One expects a 
solution path to equation (na) to be ^-Holder continuous as well, in which case 
the product Vi{zt)dXl, or the integral ^QVi{zs) dXl, cannot be given an intrinsic 
meaning since - -|- (- — l) ^0. Lyons’ deep insight was to realize that one can make 
sense of, and solve, equation ca) if one assumes one is given an enriched version 
of the driving signal X that formally consists of X together with its non-existing 
iterated integrals. The theory of regularity structures rests on the same philosophy, 
and the idea that the enriched noise should be used to give a local description of the 
unknown u, in the same way as polynomials are used to define and describe locally 
functions. 

At the very same time that Hairer built his theory, Gubinelli, Imkeller and 
Perkowski proposed in [22] another implementation of that philosophy building on 
a different notion of local description of a distribution, using paraproducts on the 
torus. The machinery of paracontrolled distributions introduced in |22] rests on a 
first order Taylor expansion of a distribution that happened to be sufficient to deal 
with the stochastic parabolic Anderson equation dLU) on the 2-dimensional torus, 
the stochastic additive Burgers equation in one space dimension |22| , the $1 equation 
on the 3-dimensional torus [laiio] and the stochastic Navier-Stokes equation with 
additive noise [38ll39]. The KPZ equation can also be dealt with using this setting 
|24] . Following Bony’s approach m , the paraproduct used in |22] is dehned in terms 
of Fourier analysis and does not allow for the treatment of equations outside the flat 
background of the torus or the Euclidean space, if one is ready to work with weighted 
functional spaces. The geometric restriction on the background was greatly relaxed 
in our previous work |1] by building paraproducts from the heat semigroup asso¬ 
ciated with the operator L in the semilinear equation. A theory of paracontrolled 
distributions can then be considered in doubling metric measure spaces where one 
has small time Gaussian estimates on the heat kernel and its ’gradient’ - see [Ij. 
This setting already offers situations where the theory of regularity structures is 
not known to be working. The stochastic parabolic Anderson model equation in a 
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2- dimensional doubling manifold was considered in [3] as an example. The first or¬ 
der ’Taylor expansion’ approach of paracontrolled calculus seems however to restrict 
a priori its range of application, compared to the theory of regularity structures, 
and it seems clear that a kind of higher order paracontrolled calculus is needed to 
extend its scope. We tackle in the present work the hrst difficulty that shows off in 
this program, which is related to the crucial use of commutator estimates between 
the heat operator and a paraproduct, which is one of the three workhorses of the 
paracontrolled calculus method, together with Schauder estimates and another con¬ 
tinuity result on some commutator. The development of a high order paracontrolled 
calculus is the object of another work [5]. 

Working in unbounded spaces with weighted functional spaces requires a careful 
treatment which was not done so far. We shall illustrate the use of our machinery 
on two examples: The parabolic Anderson model (PAM) equation (II. ID in a pos¬ 
sibly unbounded 3-dimensional Riemannian manifold, and Burgers equation with 
multiplicative noise in the 3-dimensional closed Riemannian manifold. Hairer and 
Labbe have very recently studied the (PAM) equation in from the point of view 
of regularity structures [28] - see also the work |3nj of Hairer and Pardoux. They 
had to introduce some weights w to get a control on the growth of quantities of 
interest at spatial inhnity. A non-trivial part of their work consists into tracking the 
time-behavior of their estimates, with respect to the time, which requires the use of 
time-dependent weights. For the same reason, we also need to use weighted spaces 
and working with the weights of [301 [28] happens to be convenient. Our treatment 
is however substantially easier, as we do not need to travel backwards in time such 
as required in the analysis of the reconstruction operator in the theory of regularity 
structures. As a matter of fact, our results on the (PAM) equation give an alterna¬ 
tive approach, and provide a non-trivial extension, of the results of |28| to a non-flat 
setting, with a possibly wider range of operators L than can be treated presently 
in the theory of regularity structures. As for Burgers equation with multiplicative 
noise, it provides a description of the random evolution of a velocity field on the 

3- dimensional torus, subject to a random rough multiplicative forcing, and whose 
dynamics reads 

{dt -I- L)u + [u ■ V)u = M^u, (1.3) 

where C is a 3-dimensional white noise with independent coordinates, and 

Mcu:= (CV,CV,CV), 

for the velocity field u = ^ M^. With zero noise this 3- 

dimensional Burgers system plays a very important role in the theory of PDFs 
coming from fluid mechanics, and later from condensed matter physics and statis¬ 
tical physics. R has been proposed by Burgers in the 30’s as a simplihed model 
of dynamics for Navier-Stokes equations. A change of variables, called after Cole 
and Hopf, can be used to reduce the deterministic quasilinear parabolic equation 
to the heat equation, thus allowing the derivation of exact solutions in closed form. 
Despite this fact, the study of Burgers system is still very fashionable as a bench¬ 
mark model that can be used to understand the basic features of the interaction 
between nonlinearity and dissipation. Motivated by the will to turn Burgers equa¬ 
tion into a model for turbulence, stochastic variants have been the topic of numerous 
recent works p ED Ea Ea Ea IM], where a random forcing term is added in the 
equation, mainly in one space dimension, with an additive space-time white noise 
~ that is with a space-time white noise instead of M^tt with space white noise. 
The Cole-Hopf transformation can formally be used again, and turns a solution to 
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the 1-dimensional stochastic Burgers equation with additive space-time noise to the 
heat equation with multiplicative space-time noise, with a very singular noise, such 
as detailed in [23]. A similar change of variable trick can be used for the study 
of the above multidimensional stochastic Burgers system with multiplicative noise 
we shall analyse it in Section 15.41 Also, one can consider the study of this ex¬ 
ample as a first step to understanding the dynamics of the 3-dimensional stochastic 
incompressible Navier-Stokes equation, with multiplicative noise, where the incom¬ 
pressibility brings the additional difficulty to deal with the Leray projector to keep 
the vanishing divergence property. In any case, equation seems not to have 
been studied so far, to the best of our knowledge. 

Contrary to the theory of regularity structures, whose introduction requires to 
set up a whole new algebraic-analytic setting, the analytic part of paracontrolled 
calculus is based only on classical ingredients, and its use in solving some singular 
stochastic partial differential equation involves an elementary reasoning. This ma¬ 
chinery is described in simple terms in Section [21 which serves as a baseline for the 
study of the parabolic Anderson and Burgers equations in Section O 

The geometric and functional settings in which we lay down our study are de¬ 
scribed in Section [3l In short, we work on a doubling metric measure manifold 
equipped with a Riemannian operator 

A 

i—1 

given by the finite sum of squares of vector fields. The heat semigroup of the operator 
L is assumed to have a kernel that satisfies Gaussian pointwise bounds, together with 
its iterated derivatives; precise conditions are given in the beginning of Section 13.11 
Such a setting covers a number of interesting cases. One can use the semigroup to 
construct in an intrinsic way the scale of spatial Holder spaces C°‘{M) on M and a 
scale of parabolic Holder spaces C“([0,T] x M) in which the (PAM) and Burgers 
equations will eventually be solved. Some Schauder-type regularity estimates for the 
heat semigroup, proved in Section 13.41 will be instrumental for that purpose. We 
call resolution map of the heat semigroup the map that associates to a distribution 
/ the solution to the equation (dt + L)v = /, with zero initial condition. One 
of our main contributions is the introduction of a pair of paraproducts built from 
the heat semigroup, intertwined via the resolution map, that are used to get exact 
formulas where formulas with a remainder were used previously [2211241 13]. These two 
paraproducts share the same algebraic structure and the same analytic properties, 
most importantly a cancellation property that we introduce in Section [3.21 It allows 
in particular to set the stage in a more natural function space than previously done. 
They consist in some sense of space-time paraproducts in the parabolic variable, 
whose continuity properties together with Schauder estimates allow to obtain some 
crucial estimates in L®C'"(M) spaces. 

The technical core of the paracontrolled calculus, such as defined by Gubinelli, 
Imkeller and Perkowski, is a continuity estimate for a corrector that allows to make 
sense of an a priori undefined term by compensating it by another potentially un¬ 
defined term with a simpler structure, and to separate analytic from probabilistic 
considerations. We prove in Section [4.21 that this result holds in our general setting 
as well. As a result, we are able to prove the following kind of results on the (PAM) 
in a 3-dimensional possibly unbounded measured manifold (M, d, fi) that is Ahlfors 
regular, and working with a second order differential operator L that satisfies some 
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mild assumptions stated in Section 13.11 We also study the multiplicative Burgers 
equations in a bounded ambiant space. In statements below, ^ stands for a space 
white noise on (M,/x), and := stands for its regularization via the heat 

semigroup. Full details on the mathematical objects involved in the statements will 
be given along the way. The notion of solution to the (PAM) equation (11.11) depends 
on a notion of (PAM)-enhancement C of a distribution C, e To every such 

enhancement of ( is associated a Banach space of distributions within which 

one can make sense of the equation and look for the solution to it - this is the space 
of paracontrolled distributions; see Sections [2] and O We refer to Section [331 for the 
dehnition of the weighted spaces used below, and to Section 13.11 for the statement 
of Assumptions (A) on the heat semigroup generated by L. Assumption (B) is a 
statement about the, probabilistic, renormalisation process needed to make sense of 
some ill-dehned terms; we take it for granted in the present work. It is fully spelled 
out in Section [5l and hints about the problems involved in this operation are given 
in Section O 

Theorem 1. Let (M, d,/x) be a 3-dimensional possibly unbounded metric measure man¬ 
ifold. Assume the heat semigroup satisfies Assumptions (A) and that the vector fields 
Vi are divergence-free. Let further work under assumption (B). Given , and a 

(PAM)-enhancement of a distribution ( e the parabolic Anderson model equation 

on M has a unique paracontrolled solution in T>(((). Moreover, the space white noise ^ 
has a natural (PAM)-enhancement, and there exists a sequence (-^^)o<e<i determin¬ 
istic and time-independent functions such that for every finite positive time horizon T 
and every initial data mq ^ the solution u'^ of the renormalized equation 

dtxx^ -h Lxx^ = xx= - A^), xx*'(0) = xxq 

converges in probability to the solution u e C"([0, T] x M) of the parabolic Anderson 
model equation on M associated with the natural enhancement of The result holds 
with xx; = 1 and T = oo, if /x(M) is finite. 

Let emphasize that uniqueness has to be understood as uniqueness of a solution in 
a suitable class of paracontrolled distributions, in which the problem is formulated. 
Note also that we use weighted spatial and parabolic Holder spaces to deal with 
the unbounded nature of the ambient space M. In one can typically work with 
the weights w{x,t) = and wo{x) = w{x,0) a constant - these weights were 

already used by Hairer and Labbe in [28] ; see section 13.31 Hairer and Labbe [28] are 
able to work in the range — ^ < a ^ 0, in the setting of regularity structures; we do 
not know how to deal with such a situation in our setting. Note on the other hand 
that we described in the appendix of [1] how to extend the paracontrolled calculus 
to a Sobolev setting. Together with the present work, this allows to solve the 
(PAM) equation in Sobolev spaces for a large enough finite positive exponent 

p. The above Holder setting corresponds to working with p = co. The robustness 
of our framework in terms of the operator L or the ambient geometry is useful, 
at least insofar as the tools of regularity structures have not been adapted so far 
in a non-flat setting. Moreover, as explained before, it is easier to deal with the 
time-dependent weight through the current paracontrolled approach than via the 
regularity structures theory, as done in [28] . 

As we shall see, the computations needed to handle the (PAM) equation and mul¬ 
tiplicative Burgers system involve almost the same quantities. As far as the latter is 
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concerned, we can prove the following result, under the same conditions on the op¬ 
erator L as above, in a bounded geometry. We identify in the renormalized equation 
(11.41) below a symmetric matrix d with its associated quadratic form. We can work 
in such a bounded domain with the weight w : (x,r) e^'^, for a large enough 

positive constant k. Note here that the above mentioned notion of enhancement C, 
of a distribution Q e depends on the equation under study, which is why we 

called it (PAM)-enhancement above. We denote by (B’) an assumption similar to 
assumption (B), about the renormalisation process of a number of ill-defined terms 
that appear in the analysis of the Burgers system. We work here in the 3-dimensional 
torus with valued fields u. 

Theorem 2. Assume the heat semigroup satisfies Assumptions (A) and that the vector 
fields Vi are divergence-free. Let further work under assumption (B’). Given a e (|, i), 
and a Burgers-enhancement of e (7““^, the multiplicative Burgers equation (II. 3p on 
M has a unique local in time paracontrolled solution in T’(C)- Moreover, the space 
white noise ^ has a natural Burgers-enhancement, and there exists sequences of time- 
independent and deterministic M^-valued functions (•^^)o<£;<i ^ 3)-symmetric- 

matrix-valued functions on M, such that if one denotes by the solution of 

the renormalized equation 

dtU^ Lu^ • V)u^ = — d^{u^ U^(0) = Uq (1.4) 

with initial condition UQ e then converges in probability to the solution u e 
of the multiplicative Burgers equation, locally in time. 

Details on Theorems [T] and [2] can be found in Section [5j These statements are 
two-sided, with the well-posedness of the paracontrolled version of the equations on 
the one hand, and the link between this notion of solution and the convergence of 
solutions to a renormalized regularized version of the initial equation on the other 
hand. Assumptions (B) and (B’) deal with the latter side of the stody. Note that 
after the very recent works [mills] of Bruned-Hairer-Zambotti and Chandra-Hairer 
on renormalisation within the regularity structure approach to singular PDEs, there 
is no doubt anymore that this probabilistic step should be doable in a paracontrolled 
setting as well, in some generality. 

Notations. Let us fix here some notations that will be used throughout the work. 

• Given a metric measure space (M, d, p), we shall denote its parabolic version 

by {M. wherethe parabolic space 

M := M X M 

is equipped with the parabolic metric 

p{{x,r),{y,a)) = d{x,y) -h 

and the parabolic measure n = p<S> dt. Note that for (x, t) e A4 and small 
radius r > 0, the parabolic ball ((x,r),r) has volume 

z/(B^((x,r),r)) ^ p{B{x,r)). 

We shall denote by e a generic element of the parabolic space M. 

• Given an unbounded linear operator L on we denote by 'D 2 {L) its 

domain. We give here the definition of a distribution, as it is understood in 
the present work. The definition will always be associated with the operator 
L described in Subsection 13. II below. 
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Fix a point o e M and then define a Frechet space So of test functions / 
on A4 requiring that 


Naif) ■■ = 




< 00 , 

2,du 


for all tuples a = (ai,... , 04 ) of integers; we equip So of the Frechet space 
structure associated with the following family of semi-norms Na- A distri¬ 
bution is a continuous linear functional on So', we write 5' for the set of all 
distributions. 


• Spatial Holder spaces and parabolic Space-time Holder spaces will be 
rigorously defined in Section [331 and the weights w and pa will be introduced 
in Section [331 To deal with remainder terms in some paracontrolled expan¬ 
sions, we shall use the following notation. For 7 e M and c a non negative 
integer, we shall denote by d 7 Dc an element of CZp^a ; and by (j 7 D J an element 
of 

• As a last bit of notation, we shall always denote by Kq the kernel of an 
operator Q, and write <t for an inequality that holds up to a positive mul¬ 
tiplicative constant that depends only on T. 


2 _ 

Paracontrolled calculus in a nutshell 

The theories of regularity structures and paracontrolled calculus aim at giving 
a framework for the study of a class of classically ill-posed stochastic parabolic 
partial differential equations (PDFs), insofar as they involve illicit operations on the 
objects at hand. This is typically the case in the above parabolic Anderson model 
and Burgers equations, where the products and are a priori meaningless, 
given the expected regularity properties of the solutions to the equations. So a 
regularization of the noise does not give a family of solutions to a regularized problem 
that converge in any reasonable functional space to a limit that could be defined 
as a solution to the original equation. To bypass this obstacle, both the theory 
of regularity structures and paracontrolled calculus adopt a point of view similar 
to the point of view of rough paths analysis, according to which a good notion of 
solution requires the enhancement of the notion of noise into a finite collection of 
objects/distributions, built by purely probabilistic means, and that a solution to 
the equation should locally be entirely described in terms of these objects. This 
collection of reference objects depends on the equation under study, and plays in 
the setting of regularity structures the role played by polynomials in the world of 
maps, where they provide local descriptions of a function in the form of a Taylor 
expansion. Something similar holds in paracontrolled calculus. In both approaches, 
the use of an ansatz for the solution space allows to make sense of the equation 
and get its well-posed character by deterministic fixed point methods, and provides 
as a consequence solutions that depend continuously on all the parameters in the 
problem. 

To be more concrete, let us take as an introduction to these theories the example 
of the 2 -dimensional (PAM) equation, fully studied in [20 [221 EH S]- The space 
white noise C is in that case (—l”)-Holder continuous, and the intuition suggests 
that the solution u to the (PAM) equation should be (l“)-Holder continuous, as a 
consequence of the regularizing effect of the heat semigroup. So at small time-space 







scales, u should essentially be constant, as a first approximation. This could sug¬ 
gest to try a perturbative approach in which, if one denotes by Z the solution to 
the equation {dt + A)Z = with null initial condition, one looks for a distribu¬ 
tion/function V := u — Z with better regularity than the expected regularity of u, 
so as to get a well-posed equation for v. Such an attempt is bound to fail as v needs 
to satisfy the same equation as u. The same trick invented by Da Prato-Debbusche 
in their study of the 2 -dimensional stochastic quantization equation also fails 
in the study of 3-dimensional scalar $3 equation, but a local ’version’ of this idea 
is at the heart of the theory of regularity structures, while a tilted version of that 
point of view is also the starting point of paracontrolled calculus. Both make sense, 
with different tools, of the fact that a solution should locally “look like” Z. Whereas 
’usual’ Taylor expansions are used in the theory of regularity structures to compare 
a distribution to a linear combination of some given model distributions constructed 
by purely probabilistic means, such as the a priori undehned product Z(, the para¬ 
controlled approach uses paraproducts as a means of making sense of the sentence 
“rt looks like Z at small scales”, such as given in the definition below. For readers 
unfamiliar with paraproducts, recall that any distribution / can be described as an 
infinite sum of smooth functions fi with the Fourier transform fi of fi essentially 
equal to the restriction of / on a compact annulus depending on i. A product of 
two distributions / and g can thus always be written formally as 

fa = Yi fiSj = 2 f + S + S f 

i^j-2 |i-j|<;l j^i-2 (2.1) 

= :Uf{g)+U{f,g)+Ug{f). 

The term Ilf{g) is called the paraproduct of / and g, and the term Il{f,g) is called 
the resonant term. The paraproduct is always well-dehned for / and g in Holder 
spaces, with possibly negative indices a and /3 respectively, while the resonant term 
only makes sense if a + 13 > 0. (The book [2] provides a gentle introduction to 
paraproducts and their use in the study of some classes of PDFs.) This result of 
Bony on paraproducts m already offers a setting that extends Schwartz operation 
of multiplication of a distribution by a smooth function; it is not sufficient however 
for our needs, even for the (PAM) equation in dimension 2, as u is expected there to 
be l”-Holder and C is (—l^)-Holder in that case. Needless to say, things are even 
worse in dimension 3 and for Burgers system. However, the point is that we do not 
want to multiply any two distributions but rather very special pairs of distributions. 
A reference distribution Z in some parabolic Holder space defined later, is given 
here. 

Definition. Let /3 > 0 be given. A pair of distributions {f,g) e x is said to be 

paracontrolled by Z if 

(/,ff)«:=/-n,(Z)6C“+^. 

The distribution g is called the derivative of f with respect to Z. The following 
elementary remark gives credit to this choice of name. It also partly explains why 
we shall solve the (PAM) equation in the way we do it here - using some kind of 
Cole-Hopf transform. Assume a is positive, and write (2a) for a function in 
that may change from line to line. For a pair (/, /) paracontrolled by Z, one can 
write f = e^g, for some function g in It suffices indeed to notice that Bony’s 
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decomposition gives 

e-^f = U,-z{f) + Uf{e-^) + {2a) 

= U,-z{Uf{Z)) + Uf{U_,-z{Z)) + (2a) 

= Ilg-Zf{Z) — Ilg-Zf{Z) + (2a) = (2a). 

We used in the second and third equalities two elementary results on paraproducts 
which are well-known in the classical setting, and proved below in the more general 
setting of the present work. 

The twist offered by this definition, as far as the multiplication problem of u by 
is concerned, is the following. Take for Z the solution to the equation (dt + L)Z = (, 
with null initial condition; the noise ^ is thus here (a — 2)-Holder. From purely 
analytic data, the product uC is meaningful only if a -I- (a — 2) > 0, that is a > 1. 
For a distribution [u, u') controlled by Z, with (3 = a say, the formal manipulation 

< = n„(c) + n^(w) + n(ti, c) 

= n„(c) + U^{u) + n(n„,(z), C) + n(^ 2a), c) 

=: n4C) + + C{Z, u', C) + u' U{Z, C) + n(( 2a C), 

gives a decomposition of where the first two terms are always well-defined, with 
known regularity, and where the last term makes sense provided 2a -I- (a — 2) > 0, 
that is a > |. It happens that the corrector 

C(z,u',c) :=n(n„,(z),c) -u'n(z,c) 

can be proved to define an (a-|-a-|-(a—2))-Holder distribution if a > |, although the 
resonant term n(n„/(Z), ^) is only well-defined on its own if a > 1. So we see that 
the only undefined term in the decomposition of uC, is the product u'Il{Z, (), where 
the resonant term Ii{Z,(d does not make sense so far. This is where probability 
comes into play, to show that one can define a random distribution n(Z, <^) as a 
limit in probability of renormalized quantities of the form n(Z^,<^^) — d", where 
is a regularized noise, with associated Z^ , and (f is a deterministic function, a 
constant in some cases. The convergence can be proved to hold in 
the product u' n(Z, C) eventually makes perfect sense if a -I- (2a — 2) > 0, that 
is a > |. This combination of analytic and probabilistic ingredients shows that 
one can define the product uC,^ or more properly {u,u')C, for a > |, which is 
definitely beyond the scope of Bony’s paradigm. Once the distribution ( has been 
enhanced into a pair ( := ((^,n(Z, C)) with good analytic properties, one can define 
the product {u, u') C as above for a generic distribution paracontrolled by Z, and 
reformulate a singular PDF such as the (PAM) equation in dimension 2 as a fixed 
point problem in some space of paracontrolled distribution, and solve it uniquely by 
a fixed point method. Note that the very notion of product, and hence the meaning 
of the equation, depends on the choice of enhancement of C into C- 

The above reasoning will not be sufficient, however, to deal with the (PAM) and 
multiplicative Burgers equations in dimension 3, for which a < ^, and one needs 
first to reformulate the equation differently to make it accessible to this first order 
expansion calculus. In analogy with Lyons’ rough paths theory, and in parallel with 
the logical structure of the theory of regularity structures, one may also consider 
developing a higher order paracontrolled calculus where a collection of reference 
functions {Zi, ..,Zd, with increasing regularity (for example Zi of regularity ia for 
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some a > 0), are given, and used to give some sort of Taylor expansion of a function 
/ e of the form 

for some tuple {gi, --jgk) of C“ functions with similar expansions at lower order. We 
shall develop this framework in a forthcoming work. 

3 _ 

Geometric and functional settings 

We describe in this section the geometric and functional settings in which we shall 
construct our space-time paraproducts in Section 01 and provide a number of tools. 
We shall work in a Riemannian setting under fairly general conditions; parabolic 
Holder spaces are defined Section 1,3.31 purely in terms of the semigroup generated by 
L. In Section 13.41 we prove some fundamental Schauder-type regularity estimates. 
The cancellation properties put forward in Section [3.21 are fundamental for proving 
in Section 0] some continuity results for some iterated commutators and correctors. 

3.1. Riemannian framework Our basic setting in this work will be a complete vol¬ 
ume doubling measured Riemannian manifold (M, d, /i); 
all kernels mentioned in the sequel are with respect to the fixed measure //. We are 
going to introduce in the sequel a number of tools to analyze singular partial differ¬ 
ential equations involving a parabolic operator on M+ x M 

■.= dt + L, 

with L built from first order differential operators (Vi)i=i„gg on M, defined as oper¬ 
ators Vi satisfying the Leibniz rule 

Vi{fg) = mg)+gVi{f) (3.1) 

for all functions /, g in the domain of L. Given a tuple I = (zi,..., ffc) in {1,... , ^o}^) 
we set |I| := k and 

Assumption (A) We shall assume throughout that 

• the operator L is a sectorial operator in L?‘{M), L is injective on Lp‘{M), or the 

quotient space of L?‘{M) by the space of constant functions if g, is finite, it has 
a bounded i7®-calculus on and —L generates a holomorphic semigroup 

{e~^^)t>o on L'^{M), 

• one has 

and 

V{L) c V{Vl) := {/ e ^ L\g)}, 

for some operators V) satisfies the Leibniz rule (|3.ip on 'D{L), 

• the heat semigroup is conservative, that is (e“*^)(lM) = 1 m for every t > 0, 
where 1m stands for the constant function on M - or in a weak sense that 

Li'i-Ai) = 0 , 



11 


• the semigroup has regularity estimates at any order, by which we mean that 
for every tuple I, the operators and have kernels 

Kt{x,y) satisfying the Gaussian estimate 


Kt{x,y) 


< 


. d(^,yr 


y{B{x,Vi)) 

and the following regularity estimate. For d{x, z) ^ \/t 


Kt{x,y) - Kt{z,y) 


< 


d{x, z) 


. d(x,y)-‘ 


y/t y.{B{x,y/t)) 


(3.2) 


(3.3) 


for some constants which may depend on |/|. 


Let us point out that the regularity property (13.31) for |/| = k can be obtained from 
()3.2I) with k + 1 writing the “finite-increments” formula 

\Kt{x,y) - Kt{z,y)\< d{x,z) s\XY> sup \XjKt{w,y)\ 

j wg(x,z) 

where {x,y) stands for a geodesic joining x to z and of length d{x,z), and (Xj) 
stands for a local frame field near (x,y), and it acts here as a first order differential 
on the first component of K. As a matter of fact, it suffices for the present work to 
assume that the semigroup has regularity estimates of large enough order. Observe 
that under Assumption ??, the semigroup may be defined as acting on the 
distributions. This can be rigorously done by duality, since for every integer ^ 0, 
[L *)^has a kernel satisfying pointwise Gaussian estimates and so is acting on 
the test functions So - we refer to | 12 ] for more details. One can keep in mind the 
following two examples. 


(a) Euclidean domains. In the particular case of the Euclidean space, all of 
the current work can be reformulated in terms of Fourier transform rather 
than in terms of the heat semigroup; which may make some reasoning a bit 
more familiar but does not really simplify anything. The case of a bounded 
domain with its Laplacian associated with Neumann boundary conditions fits 
our framework if the boundary is sufficiently regular. We may also consider 
other kind of second order operator, like L = — div(AV) for some smooth 
enough matrix-valued map satisfying the ellipticity (or accretivity) condition. 


(b) Riemannian manifolds. Smooth closed manifolds equipped with an oper¬ 
ator L of Hormander type as above, with Vi smooth, all satisfy assumptions 
(A). Here is a simple setting within which one can deal with unbounded 
spaces. Assume M is a parallelizable d-dimensional manifold with a smooth 
global frame field V = (Hi,...,Vd). One endows M with a Riemannian 
structure by turning V into orthonormal frames. The above assumption on 
the heat kernel holds true if M has bounded geometry, that is if 

(i) the curvature tensor and all its covariant derivatives are bounded in the 
frame field V, 

(ii) Ricci curvature is bounded from below, 

(iii) and M has a positive injectivity radius; 

see for instance m or m- One can include the Laplace operator in this 
setting by working with its canonical lift to the orthonormal frame bun¬ 
dle, given by 5 where the Hi are the canonical 
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horizontal vector fields of the Levi-Civita connection, and the Vjk are the 
canonical vertical vector fields on the orthonormal frame bundle, inherited 
from its 5'0(M'^)-principal bundle structure. The bundle OM is parallelizable 
and satisfies assumptions (A) if the base Riemannian manifold M satisfies 
the conditions (i-iii). This example shows that the assumption that M is 
parallelizable is essentially done here for convenience. 


3.2. Approximation operators and cancellation We introduce in this section a notion 
property of approximation operators that will 

be the building blocks for the defini¬ 
tion and study of the paraproducts, commutators and correctors, used in our analysis 
of singular PDEs. Some of them enjoy some kind of orthogonality, or cancellation, 
property quantified by condition (l3.1Up below. This property is to be thought of as 
a quantitative replacement for the property of frequency localisation of Fourier mul¬ 
tipliers in the Littlewood-Paley decomposition heavily used in the classical Fourier 
definition of paraproducts; see [2]. Indicators of annulii will somehow be replaced 
in our setting by continuous functions of order 1 on such annulii, with exponential 
decrease at 0 and oo. Note that we shall be working in a parabolic setting with 
mixed cancellation effects in time and space. 

All computations below make sense for a choice of large enough integers h, 
that will definitely be fixed at the end of Section 14.11 to ensure some continuity 
properties for some useful operators. Recall that generic elements of the parabolic 
space A4 = M X M are denoted by e = (x, r) or e' = {y, a), and that t stands for a 
scaling parameter. The following parabolic Gaussian-like kernels {Gt)o<tsii will be 
used as reference kernels in this work. For 0 < t ^ 1 and a ^ r, if d{x, y) ^ 1, set 


Gt{{x,T),{y,a)) : = 




1 + 


p{{x,T),{y,a)) 


2\ -^1 


otherwise we set 
Gt{{x,T),{y,a)) := 
1 




1 + 


r — cj 


-ei 


d{x,yf 


2 \ -^1 


exp —c 


d{x,yf 


for d{x,y) ^ 1, and = 0 if r ^ cj. We do not emphasize the dependence of G on 
the positive constant c in the notation for the ’Gaussian’ kernel, and we shall allow 
ourselves to abuse notations and write Gt for two functions corresponding to two 
different values of that constant. This will in particular be the case in the proof of 
LemmalU We have for instance, for two scaling parameters s,t e (0,1), the estimate 

r Gt (e, e') Gs (e', e") i/{de') < Gt+s (e, e"). (3.4) 

JM 

(Indeed, the space variables and the time variables are separated in the kernel Gt- 
Then both in space and time variables, the previous inequality comes from classical 
estimates for convolution of functions with fast decay at infinity, such as done in [H 
Lemma A.5] for example.) This somewhat unnatural definition of a Gaussian-like 
kernel is justified by the fact that we shall mainly be interested in local regularity 
matters; the definition of Gt in the domain |d(x,y) ^ 1} is only technical and will 
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allow us to obtain global estimates with weights. Presently, note that a large enough 
choice of constant l\ ensures that we have 

sup sup Gt{G,e') ^{de') < CO, (3.5) 

ls(0,l] eeM JM 

so any linear operator on a function space over M, with a kernel pointwisely bounded 
by some Gt is bounded in for every p e [1, oo]. 


Definition. We shall denote throughout by G the set of families ('Pt)o<isci of linear 
operators on M with kernels pointwisely bounded by 


K-pt{e,e') 


< Gtie,e'). 


The letter G is chosen for ’Gaussian’. A last bit of notation is needed before we 
introduce the cancellation property for a family of operators in a parabolic setting. 
Given a real-valued integrable function m on M, define its rescaled version as 


m-t(-) 



the family (mt)o<t^i is uniformly bounded in L^(M). We also define the “convolu- 
tion” operator m* associated with m via the formula 

^00 

:= m{T - a)f[a)d(T. 

Jo 


Note that if m has support in M+, then the operator m* has a kernel supported 
on the same set {{a,T);a ^ r} as our Gaussian-like kernel. Moreover, we let the 
reader check that if mi, m 2 are two L^-functions with m 2 supported on [0, 00 ), with 
convolution mi * m 2 , then we have 


(mi * m 2 )* = m^ o m^. 

Given an integer b ^ 1, we define a special family of operators on L?‘[M) setting 
76 := (5 — 1)! and 


■= lb 


^-tL 


and 




with = Id, so is an operator of the form pb{tL)e~^^, for some polynomial 
Pb of degree (5 — 1), with value 1 in 0. Under Assumption ??, the operators Pj: 
and both satisfy the Gaussian regularity estimates (|3.2I1 at any order 


( 6 ) 


K |/| (x,y) 

V 

K \i\ (x,y) 

t^ViR 


t^RVi 


< 


. d{x,y)^ 




(3.6) 


with R standing here for or Q 


lb) 


The parameters 5 and li will be chosen large enough, and fixed throughout the 
paper. See Proposition [15] and the remark after Proposition [16] for the precise choice 
of 5 and li. 

Definition. Let an integer a e |0, 25] be given. The following collection of families of 
operators is called the standard collection of operators with cancellation of order 

a, denoted by StGC“. It is made up of all the space-time operators 
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where k is an integer with 2k + \J\ ^ o, and c e [1, 6], and m is any smooth function 
supported on [ 5 , 2 ] such that 

J T'‘m{T) dr = 0, (3.7) 

for all 0 ^ i ^ k — 1, with the first b derivatives bounded by 1. These operators are 
uniformly bounded in for every p e [l,oo], as functions of the scaling parameter 

t. So a standard collection of operators Q can be seen as a bounded map Q '■ t ^ Qt 
from (0,1] to the set B{U’) of bounded linear operators on We also set 

StGCP’^^] := [j StGC“. 

0^a^2b 

The cancellation effect of such operators is quantified in Proposition [3] below; 

a — I J\ — 2k 

note here that it makes sense at an intuitive level to say that L 2 encodes 
cancellation in the space-variable of order a — | J| —2k, that Vj encodes a cancellation 
in space of order \J\ and that the moment condition (13.7|) encodes a cancellation 
property in the time-variable of order k for the convolution operator m*. Since 
we are in the parabolic scaling, a cancellation of order k in time corresponds to a 

a~|J|— 2 fc ( 

cancellation of order 2k in space, so that VjL 2 0 rn* has a space-time 

cancellation property of order a. We invite the reader to check that each operator 

/ 1^1 \ . s Q—I J| —2fc 

(^t 2 Vjj [tL) 2 ^ (S) Tn* in the standard collection has a kernel pointwisely 

bounded from above by some Qt- This justifies the choice of name StGC“ for this 
space, where St stands for ’standard’, G for ’Gaussian’ and C for ’cancellation’. 
The paracontrolled analysis, that we are going to explain, is based on these specific 
operators. We emphasize that because of the Gaussian kernel Gt and the function 
m, all of these operators have a support in time included in 

{(r,cr), T ^ fj}. 

In particular, that means that we never travel backwards in time through these 
operators. This fact will be very important, to deal further with the weight w, which 
will depend on time. We give one more definition before stating the cancellation 
property. 

Definition. Given an operator Q := Vj 4>{L), with |/| ^ 1, defined by functional calculus 
from some appropriate function cj), we write Q* for the formal dual operator 

Q* := ^{L)Vj. 

For 1 = 0, and Q = (l){L), we set Q' := Q. For an operator Q as above we set 

[Q0m*y := Q* 0 m*. 


Note that the above definition is not related to any classical notion of duality, and 
let us emphasize that we do not assume that L is self-adjoint in L^(^). This notation 
is only used to indicate that an operator Q, resp. Q*, can be composed on the right, 
resp. on the left, by another operator 'G{L), for a suitable function i/j, due to the 
functional calculus on L. In the setting of analysis on a finite dimensional torus, the 
operators are given in Fourier coordinates A, as the multiplication operators 
by as this function is almost localized in an annulus |A| ~ t~^, the 


operators and are almost orthogonal if | is either very small or very big. 
This is encoded in the elementary estimate 

ts 




< 


[s + ty) \/s + t)) 


exp —c 


d^{x,y) 
t -\- s 


(3.8) 
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The frequency analysis of the operators is not very relevant in the non-homogeneous 
parabolic space A4. We keep however from the preceeding analysis the idea that 
relation (j3.8p encodes some kind of orthogonality, or cancellation effect. 


Proposition 3. Consider Q} G StGC"^ and G StGC“^ two standard collections with 
cancellation, and set a := min(ai,a 2 ). Then for every s,t e (0,1], the composition 
Q\ ° 2i* ^ kernel pointwisely bounded by 


K. 




.(e,eO 


< 


((s + tp) 


2 




(3.9) 


Proof - Given 

QI = ~0ml* and Qf* = ~0 

a standard operator and the dual of another, we have 

Qj O Qr = 0 {ml * m?)*. 

Assume, without loss of generality, that 0 < s ^ t. Then the kernel of the 


time-convolution operator * m[^^ is given by 




i*^2(t - ct) = ("T”) 


X — a\ dX 
st 


Since m^ has vanishing ki hrst moments, we can perform ki integration by parts 
and obtain that 

I f T — X\ r, f X — a\ dX 


K„ 


' 1 ifcrr)^ 


{t^cf) 


< 


S\ ^1 

l) 


I 


d-kl 


m 




t 


st 


where we slightly abuse notations and write d ^^m^ for the kf^ primitive of m^ 
null at 0. Then we get 




|t - A| 


-£i+2 


1 + 


X — a 


-ei+2 


st 


<^S^kF 


1 + 


|r — a\ 
s + t 


-h 


{s + t) 


-1 


_ _ ... . _ _ ~J 1 ~ 2 fc-|-l-a^ —_ 

In the space variable, the kernel of Vj^L 2 p^ is bounded 

above by 


{s -I- 1) 


— cij — ci‘2 -i~'2k‘2 


m B{x, \/s + t 


exp — c 


s + t 


as a consequence of the property (13.6p . Altogether, this gives 

/S\kl -2fel a2~2fc9 . -til +2^1 -a2 + 2fc2 , . 

^SioQ2*(e,e) ^(7) S 2 t 2 (s + t) 2 gt+s{e,e) 


< 


< 


s\ — 
t 


Ot+s (®T ) 
(^) " St+s{e,e'), 


where we used that s ^ t and a ^ ai. 


\> 
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Definition. Let 0 ^ o ^ 26 be an integer. We define the subset GC“ of G of families 
of operators with the cancellation property of order a as the set of elements Q of 
G with the following cancellation property. For every 0 < ^ 1 and every standard 

family S s StGG“ , with a' e |a,26]], the operator Qt o S* has a kernel pointwisely 
bounded by 


^QtoS’{e, 




2 


(3.10) 


Here are a few examples. Consider a smooth function m with compact support 
in [2“^, 2], an integer c ^ 1, and a tuple I of indices. 

• The families and (Vj0 belong to GG“ if 

V /0<ts:i V /0<fsci 

l-^l ^ a; 

• If 5 T^m{T) dr = 0 for all integer A; = 0,..., a—1, then we can see by integration 
by parts along the time-variable that (p/' 0 ™'*)o<r<i ^ 

• If ^T^m{T) dr = 0 for all integer k = 0, ...,02 with ai -I- 02 = o, then the 

families 0mfl and 0mll , where |/| ^ ai, both 

V /o<r^i V /o<r!£i 

belong to GC“. 


We see on these examples that cancellation in the parabolic setting can encode 
some cancellations in the space variable, the time-variable or both at a time. 


We introduced above the operators and acting on the base manifold M. 
We end this section by introducing their parabolic counterpart. Choose arbitrarily 
a smooth real-valued function </? on M, with support in [^, 2 ], unit integral and such 
that for every integer /c = 1 ,..., 6 , we have 


J T^Lpij) dr = 0 . 


Set 

Pf) := P/') 0 yp* and Q? := ^ , 

Denote by M,- the multiplication operator in M by r. An easy computation yields 
that 

sf 0<^* + pf)0V^i 

where ■= (For an extension of the present theory to the setting of 

Sobolev spaces, such as done in the appendix B of [1], it would be convenient to work 
with (p * (p rather than ip.) Note that, from its very definition, a parabolic operator 
belongs at least to GC^, for 6^2. Remark that if C is a time-independent 

distribution then Q.f'\ = Note also that due to the normalization of ip, then 

for every / 6 LP(M) supported on [ 0 , 00 ) then 


So, the operators Vt tend to the identity as t goes to 0, on the set of functions 
/ e L^{M) with time-support included in [ 0 , 00 ), whenever p e [l,oo), and on the 
set of functions / e C^{M) with time-support included in [0,oo). The following 
Calderon reproducing formula follows as a consequence. For every continuous 
function / e L°°(AI) with time-support in [ 0 ,oo), we have 

/= f 

Jo ^ 


(3.11) 
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This formula will play a fundamental role for us. Noting that the measure ^ gives 
unit mass to intervals of the form 2“*], and considering the operator as 

a kind of multiplier roughly localized at frequencies of size Calderon’s formula 
appears as nothing else than a continuous time analogue of the Paley-Littlewood 
decomposition of /, with ^ in the role of the counting measure. 


3.3. Parabolic Holder spaces We define in this section space and space-time weighted 

Holder spaces, with possibly negative regularity index, 
and give a few basic facts about them. The setting of weighted function spaces is 
needed for the applications to the parabolic Anderson model and multiplicative 
Burgers equations on unbounded domains studied in Section [5l The weights we use 
were first introduced in [28] . 

Let us start recalling the following well-known facts about Holder spaces on M, 
and single out a good class of weights on M. A function w : M ^ [IjOo) will be 
called a spatial weight if one can associate to any positive constant ci a positive 
constant C 2 such that one has 

^C2w{y), (3.12) 


for all X, y in M. Given 0 < a ^ 1, the classical metric Holder space is defined as 
the set of real-valued functions / on M with finite Lf"-norm, defined by the formula 


\f\\HS ■= Ik V 




\f{x) - f{y)\ 

sup ' , ^ 

0<d{x,y)^l w{x)d{x,y)^ 


< 00 . 


Distributions on M were defined in [3] using a very similar definition as in the end 
of Section [H where their parabolic counterpart is defined. 


Definition. For a e (—3,3) and w a spatial weight, define C" := C^{M) as the set of 
distributions on M with finite C^-norm, defined by the formula 


lla“ := w k 


-F sup t 2 
L^{M) o<t^l 


w-^Q^f 


L^{M) 


and equip that space with the induced norm. The latter does not depend on the integer 
a > and one can prove that the two spaces and C" coincide and have equivalent 
norms when 0 < a < 1 - see 0. 


These notions have parabolic counterparts which we now introduce. A space- 
time weight is a function lo : A4 ^ [l^co) with uj{x, •) non-decreasing function of 
time, for every x e M, and such that there exists two constants ci and C 2 with 

Uj(x, t) ^ C 2 Cj(y, r), (3.13) 


for all pairs of points of A4 of the form ((x, r), (y, r)). The function Wr '■= w(-,t) 
is in particular a spatial weight for every time r. For 0 < a ^ 1 and a space-time 
weight cj, the metric parabolic Holder space 7^2 = ^S(-^) is defined as the set of 
all functions on Ad with finite defined by the formula 


I/Iks := Ik V 


L^{M) 


sup 

0<p({x,T),{y,a)) SCI; 


\f{x,T) - f{y,a)\ 

lo{x,t) p{{x,T),{y,a)y 


As in the above spatial setting, one can recast this definition in a functional setting, 
using the parabolic standard operators. This requires the use of the following ele¬ 
mentary result. Recall that the kernels Gt depend implicitly on a constant c that 
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may take different values with no further mention of it. We make this little abuse 
of notation in the proof of this statement. 

Lemma 4. Let A be a linear operator on A4 with a kernel Ka pointwisely bounded by 
a Gaussian kernel Qt, for some t e (0, 1]. Then for every space-time weight uj, we have 


Proof — Indeed, for every (x, r) e Ad we have 
1 


Lo{x, r) 


JM 

f Gt{ix,T),{y,a)) 
JM 

f Ot((x,T),(y,cr)) 
JM 

< ||a;“^/|| , 

nn * 


< 


< 


oj(x,t) oj(y,a) 
u;(x,(t) CAj(y,a) 


v[dyda) 

i'{dyda) 


uj{y,a) 


v{dyda) 


where 

• we used in the second inequality the fact that the function iu(x, •) of time 
is non-decreasing, and Qt is null if cr ^ r, 

• the implicit constant in Qt was changed in the right hand side of the third 
inequality, and we used the growth condition (|3.13p on w as a function of 
its first argument here, 

• we used the uniform bound p.Sp on a Gaussian integral in the last line. 

[> 


Recall that distributions were introduced in the end of Section [TJ 

Definition. For a e (—3,3) and a space-time weight uj, we define the parabolic Holder 
space := C°(A4) as the set of distributions with finite C“-norm, defined by 

ll/llcg := sup ||W^Qi(/)||^oo(_m) + 

SeStGC*: ^ geStGC*: 0<1^1 ^ ^ 

0^k^2b |Q:|<fc^26 

equipped with the induced norm. 

The restriction a e (—3,3) is irrelevant and will be sufficient for our purpose in 
this work; taking b large enough we can allow regularity of as large an order as we 
want. Building on Calderon’s formula (13.lip , one can prove as in [1] that the two 
spaces 7^2 and C“ coincide and have equivalent norms, when 0 < a < 1. 

Proposition 5. For a e (0,1) and every space-time weight uj, the two spaces and 
C" coincide and have equivalent norms. 

Proof — We hrst check that is continuously embedded into C“. So fix a function 
/ e T-L ^, then by Lemma H] we easily deduce that 

sup ||W^Ql(/)||^o.(^) < ||w■VLao(^)• 

QeStGC^ 

0^k^2b 
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For the high frequency part, we consider t e (0,1] and Q e StGC^ with a < k ^ 
2b. Then Qt has at least a cancellation of order 1, hence 

Qt{f){e) = Qt{f - f{e)){e) 

= J Kq, (e, e') {f{e') - /(e)) i^{de'). 

Due to the kernel support of Qt, the integrated quantity is non-vanishing (and 
so relevant) only for t ^ a, with e = {x,t) and e' = {y,(7). If p{e,e') ^ 1, then 
by definition 

|/(e0 - /(e)| ^ a;(e)p(e,e')“l|/||ws 

and if p{e',e) ^ 1, then by the property of the weight we have 
|/(eO-/(e)| ^ {u{e)+uj{e'))\\u-^f\\^^^j^y 

Hence 

\Qtif){e)\ <a;(e)|J Qt{e,e')p{e,e'Y iy[de') 

< u{e)t^\\f\\u^, 

uniformly in e e Ad and t e (0,1); this concludes the proof of the continuous 
embedding of 7^2 into C“. 


To prove the converse embedding, let us start by fixing a function / e C“. The 
low frequency part of / is easily bounded, using Calderon’s reproducing formula 


11 ^ A\l<^{M) ~ 


^ !I/IIc 2 > 


L^(M) Jo 


r 


cu-IqS'V 


dt 

L”(M) t 


since a > 0. Now fix e = (x, r) and e' = {y, a) in Ad, with p := p{e, e') ^ 1 and 
T ^ (T. We again decompose 


f = V 




For t < p'^, we have 


QrV(e) <t 2 u{e)\\f\\cc 


and 


Qi^V(eO < t2w(e')||/||cs ^t2ui{e)\\f\\c^ 


where we used that the weight is increasing in time and then that d{x, y) ^ p 
with the property of the weight. So we may integrate over t < p^ and we have 

QrV(e) - QrV(e') j < f w(e)||/||c2 


< 


/3“w(e)||/||cs- 
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For the low frequency parts, with ^ t ^ 1 or 'P^\ we use that 


Q^^fix.T) - Qf>f{x,a) 


( 1 ). 


< T — (T h 


sup 

,?e(cr,r) 


drdf' f{x,Ci) 


sup 

,?e(o-,T) 




< pu){x,T)t 2 \\f\W 


where we used that p ^ 1 with the fact that the two collections of operators 
and are of type StGC^, that is have cancellation of 

order at least 1, and that the weight is non-decreasing in time. Similarly we can 
estimate the variation in space with the assumed finite-increment representation 
(lOD . where one considers a local frame field {Xj) in a neighbourhood of a 
geodesic (x, y) from x to y. This gives 


f{x,a)-Qf’f{y,a) <d{x,y) sup XjQf‘f{z,c;) 

ze{x,y) -J 


l(l) 


FI) 


< pu;{x,T)t 2 ll/llco. 


So we get 


L - ®™/(e')| j S P (J (";■ j) ^{e)\\fk 


< p‘^u}{e) ll/llcs, 

because a < 1. A similar estimate for ends the proof of continuous embed¬ 
ding of into 


The next proposition introduces an intermediate space whose unweighted version 
was first introduced in the setting of paracontrolled calculus in [22] , and used in [3] . 

To fix notations, and given a space-time weight a;, we denote by [Ct = 

{L^Ct )(a;) the set of parabolic distributions such that 


sup ||/(x,-)lLf 

XGAI uj{x,-) 


(R+) 


< 00 . 


Also 


(cj) stands for the set of parabolic distributions such that 


sup < 00 . 


Proposition 6. Given a e (0,2) and a space-time weight cu, set 

C:= (c|L®)(a;)n(L“C“)(a;). 

Then T" is continuously embedded into C". Furthermore, if a e (0,1), the spaces 
and T-L^ are equal, with equivalent norms. 


Proof — We first check that T" is continuously embedded into C", and fix for that 
purpose a function / e T". As done in jU Proposition 2 . 12 ], we know that for 
all integers k,j with k + ^ ^ and every space function g e C^{M), we have 


tWj{tLfe-^^g 




< t2 


g\\c°‘iM), 
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for any subset of indices J with \J\ = j. So consider a generic standard fam¬ 
ily 0 m*) in StGC“, with 3 ^ a ^ 6, and a smooth 

function m with vanishing first k moments. If A; = 0 we have seen that we have 




L'»(M) 


^t2\\f{T)\ 


C° 


for every r, so 






ct ,, 


L®C“(a;) 


since m* is a L®(M)-bounded operator as a convolution with an L^-normalized 
function. 

If A: = 1 (or A: ^ 1), the same reasoning shows that we have 


io{x,-) ^m*{f){x,-) 


L®(K+) 




"C 


i^r) 


for every x e M, since ^ e (0,1), and m encodes a cancellation at order 1 in 
time as it has a vanishing hrst moment. Hence 




L^(M) 






which concludes the proof of the embedding ^ C“. The remainder of the 
statement is elementary since C“ = is embedded in £^". 

[> 


Before turning to the dehnition of an intertwined pair of parabolic paraproducts, 
we close this section with two other useful continuity properties involving the Holder 
spaces C^. 

Proposition 7. Given a e (0,1), a space-time weight uj, some integer a ^ 0 and a 
standard family P e StGC“, there exists a constant c depending only on the weight w, 
such that 


^ ('Ptf)(e) - (Psf)(e') < (s+ t + p(e,e')^)^ 


uniformly in s,t e (0,1] and e = (x, r) and e' = (y,a) e M, with t ^ a. 

Proof — We explain in detail the most difficult case corresponding to P e StGC®, so 
V encodes a priori no cancellation. Then Vt takes the form 

Vt = 0 m* 


for some integer c ^ 1 and some smooth function m. There is no loss of gen¬ 
erality in assuming that ^m^r) dr is equal to 1, as V is actually an element of 
StGC^ if m has zero mean - this case is treated at the end of the proof. 


In this setting, since / is bounded and continuous, we have the pointwise identity 

/ = 
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(i) Consider first the case where p{e,e') ^ 1„ with e = (x,r) and e' = (y, cr). 
Decompose 

w(e)"^ {Vtf){e) - {Vsf){e') 

^ {'Ptf){e) - /(e) 

+ a;(e)“^ {Vtf){e) - f{e') 

< u{e)-^\f{e) - f{e')\ + ||a;"^ {Vtf - f) [Vsf - f) Iloo^mY 

We have 

^ P(e,e')"||/||ws ^ P(e, e')"||/||cs- 

For the two other terms, we use that 


L^{M) 


< 


f 


w ^uduVuf 


du 

L^iM) U ’ 


and note that 

uduVu = ®mu + Pjf'^ ® ku 

with A:(r) = dr(Tm(r)), is actually the sum of two terms in StGC^^ since it is 
clear for the hrst one and the function k has a vanishing hrst moment. It follows 
by dehnition of the Holder spaces with a <1, that we have 

- f)\\Lao^M) ~ v) ll/llcs ^i^ll/llcs- 

A similar estimate holds by replacing t by s, which then concludes the proof in 
this case. 

(ii) In the case where p{e^ e') ^ 1, we do not use the difference and use condition 
(|3.13p on the weight oj to write 

and obtain as a consequence the estimate 

^(e)"^ {Ptf){e) - {Vsf){e') ^ w(e)"^ {Ptf){e) + u}{e)~^ {Psf){e) 


< 






for some positive constant c. Since we know by Lemma 0] that Vt and Vs are 
bounded in L“(w), we deduce that 


a^(e) ^ {Ptf){e) - {Vsf){e') 


^ e' 


cd{x,y)^^-l jj 




^ „cd(x,y) II jI 

since C" L“, given that a > 0. The expected estimate follows from that 
point. 

• In the easier situation where V e StGC“ for some integer a ^ 1, we can perform 
the same reasoning and use in addition the fact that 

lim VtU) = 0, 

which makes the case easier since we do not have to deal with the first term 
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[> 


With an analogous reasoning (indeed simpler) we may prove the following. 
Proposition 8. Given a e (—3,0), a space-time weight uj and a standard family P e 
StGC°, one has 

\\'Ptf\\L‘^{M) ~ ll/llcs; 

uniformly in t e (0,1]. 

Proof — The proof follows the same idea as the the proof of Proposition [71 Indeed, 
we use the fact that since P is a standard family then 

d s 

Vtf= {-sdsPs)f- + Pif. 

Jt s 

The key point is that {—sdsPs)s can be split into a finite sum of families of 
StGC^^, which allows us to conclude as previously. 

[> 


3.4. Schauder estimates We provide in this subsection a Schauder estimate for the 

heat semigroup in the scale of weighted parabolic Holder 
spaces. This quantitative regularization effect of the heat semigroup will be instru¬ 
mental in the proof of the well-posedness of the parabolic Anderson model (PAM) 
and multiplicative Burgers equations studied in Section O Define here formally the 
linear resolution operator for the heat equation by the formula 

^-\v)r:=f (3.14) 

Jo 

We fix in this section a finite positive time horizon T and consider the space 

A4t ■= M X [0,T], 

equipped with its parabolic structure. Denote by L® the corresponding function 
space over [0,T]. We first state a Schauder estimate that was more or less proved 
in the unweighted case in [221 0 ] - see Lemma A.9 in |22] and Proposition 3.10 in 

In- 

Proposition 9. Given /3 e (—2,0) and a space-time weight to, we have 

\\^-\v)\\c,+2 <T 

We shall actually prove a rehnement of this continuity estimate in the specific 
case where to has a special structure motivated by the study of the (PAM) and 
multiplicative Burgers equations done in Section (5) These special weights were hrst 
introduced by Hairer and Labbe in their study of the (PAM) equation in and , 
via regularity structures [271[28]. Let o = Oj-ef be the reference point in M hxed and 
used in the dehnition of So at the end of Section [TJ and set 

Pa(a;) := (l-h d(oref,x))“, n7(3:,T) := (3.15) 

for 0 < a < 1 and a positive constant k. (The introduction of an extra exponential 
factor in our space-time weight w will allow us to get around an iterative step 
in the forthcoming application of the hxed point theorem used to solve the (PAM) 
and multiplicative Burgers equations, as done in [ZaEH].) For r ^ 0, we use the 
notation 


ZUr : X 6 M t— > w{x, t) 
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for the spatial weight. The space-time weight w satisfies condition (13.131) on [0, T] x 
M, uniformly with respect to k > 0. The above special weights satisfy in addition 
the following crucial property, already used in [27112H] • We have 

Pa{x)w{x,cr) < k~^{t — w{x,t), (3.16) 

for every non-negative real number e small enough, uniformly with respect to x 6 
M, K > 0 and 0 < cr < r ^ T. The next improved Schauder-type continuity estimate 
shows how one can use the above inequality for the specihc weights to compensate 
a gain on the weight by a loss of regularity. 

Proposition 10. Given /3 6 M, a 6 (0,1) and e 6 [0,1) small enough such that a+e < 1, 
we have the continuity estimate 

11 -^ ~ ^ M{L^C^){zupay 

Moreover if —2 -I- 2(a -I- e) < /3 < 0, then 


Proof — Let us hrst check the regularity in space. So consider an integer c ^ -I-1 

and a parameter r e (0,1]. Then for every hxed time r e [0, T] we have 

Jo 

By using the specihc property (|3.16l) of the weights pa and w, one has 


L^(M) 


< K 


r + T — a 

C 








r + T — a 


{r + T - a )2 {t - a) “ 


So by integrating and using that c is taken large enough, we see that 


< K~ 




This holds uniformly in r e (0,1] and r e [0,T] and so one concludes the proof 
of the hrst statement with the global inequality 


^{v). 


< K 




^da\ II u 


L'»(M) ' (Jo ' ' J " 

~ (L”Cf) (paro)' 

For the second statement, we note that for 0^iT<r^rwe have 
^-\v)r - - Id)^-i(u), + J" dr 

= r —+ re-("-")Vhr. 

Jo r 
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W. 


We have by the previous estimate 
dr 
r 


r 

Jo 


< K 


L^{M) 


i 


" \^+l-a-e±_ 


v)cr \\^l3+2~2a-2e 
II O-ci- 




(L”Cf ) {m) 


where we used that Wr ^ rua for cr ^ r. Moreover, since j3 is negative, we also 
have 




■‘f 




dr 


< 


£(£ 

r 




T - r 




P 'Vr 

1 


ds 

L'^(M) S 




{vr) 




dr 


(T-r) “ r s2—+ (T-r)’ 
Jt- r S 

where we used (I3.16P and ^ + 1 — a — e>0. 


„-L 


(vr) 


dr 


> 


The following result comes as a consequence of the proof, combined with Lemma 
m we single it out here for future reference. 

Lemma 11. Let ^ be a linear operator on A4 with a kernel pointwisely bounded by Gt 
for some t e (0,1]. Then for every a + e e (0,1), we have 

\\MLSpJM)^L<SiM) ^ 

Schauder estimates can also be extended to spaces of positive regularity. 

Proposition 12. Given P e (0,2), a e (0,1) and e e [0,1) small enough such that 
a + e < 1, we have the continuity estimate 

\J^~^{v)\„p+2-2a-2e < K“^|uL/9 . 

Proof - This follows from Proposition fTOl For v e Cmpa ^ {L'^Cx){'cupa), it 
is known that Lv e (^L^Cx~^) (^Pa), to which Proposition [10] can be ap¬ 
plied since /3 — 2 < 0. Now use that and L commute to deduce that 

L(Jf~^v) e hence Jf~^v e (ro). On the other hand, 

dt(^~^v) = V — from which follows that dt^~^v e and con¬ 
sequently e 

[> 

The constraint /3 < 2 is not relevant. Indeed, by iteration the previous Schauder 
estimates can be proved for an arbitrary exponent /3 > 0. 


4 _ 

Time-space paraproducts 

We introduce in this section the machinery of paraproducts which we shall use 
in our analysis of the singular PDFs of Anderson p.ip and Burgers (|1.3I) . In the 
classical setting of analysis on the torus, the elementary definition of a paraproduct 
given in Section [2] in terms of Fourier analysis should make convincing, for those 
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who are not familiar with this tool, the fact that nj((5r) is a kind of ’’modulation” 
of g, insofar as each mode gj of g, in its Paley-Littlewood decomposition, is mod¬ 
ulated by a signal which oscillates at frequencies much smaller - the finite sum 
fi- makes sense to talk of a distribution/function of the form Ilf{g) 

as a distribution/function that ’’locally looks like” g. This is exactly how we shall use 
paraproducts, as a tool that can be used to provide some kind of Taylor expansion 
of a distribution/function, in terms of some other ’model’ distributions/functions. 
This will be used crucially to bypass the ill-posed character of some operations in¬ 
volved in the (PAM) and Burgers equations, along the line of what was written in 
Section [2l 

Working in a geometric setting where Fourier analysis does not make sense, we 
shall define our paraproduct entirely in terms of the semigroup generated by the 
operator ^ = dt + L on the parabolic space. The dehnition of a paraproduct comes 
together with the dehnition of a resonant operator n(-, •), tailor-made to provide the 
decomposition 

fg = Uf{g)+U{f,g)+Ug{f) 

of the product operation, and with Ilf{g) and n(/, g) with good continuity properties 
in terms of / and g in the scales of Holder spaces. Such a construction was already 
done in our previous work [3], where the generic form of the operator L, given by 
its Hrst order carre du champ operator, imposed some restrictions on the range of 
the method and allowed only a hrst order machinery to be set up. The fact that we 
work here with an operator L in Hormander form will allow us to set up a higher 
order expansion setting. We will use this for the description of the space in which 
to make sense of the two singular PDFs we want to analyse. However, this a priori 
useful setting is in direct conhict with one of the main technical tools introduced by 
Gubinelli, Imkeller and Perkowski in their seminal work |22j . 

The case is easier to explain on the example of the (PAM) equation. A solution 
to that equation is formally given as a hxed point of the map 

$ : u e~'^UQ + 

for which we shall need tt to be a priori controlled by Z := to make sense 

of the product uC, - more will actually be required, but let us stick to this simplihed 
picture here; so the map $ will eventually be defined on a space of distributions 
controlled by Z, such as dehned in Section [21 where it will be shown to be a con¬ 
traction. At a heuristic level, for a distribution (tt, u') controlled by Z, the product 
uC, will be given by a formula of the form 

nc = n„(C)+ (•••)• 

To analyse the term and recalling that Z := it is thus very 

tempting to write 

if-i(n„(c)) = u^z) + [^-\Uu]{C) + (•••) 

and work with the commutator This is what was done in [23 m to 

study the 2-dimensional (PAM) equation on the torus and more general settings; 
and it somehow leads to a non-natural choice of function space for the remainder 
of a paracontrolled distribution in a space-time setting. Unfortunately, we have 
little information on this commutator, except from the fact that it is a regularizing 
operator with a quantifiable regularizing effect - it was hrst proved in [22] in their 
Fourier setting. This sole information happens to be insufficient to push the analysis 
of the (PAM) or Burgers equations far enough in a 3-dimensional setting. As a way 
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out of this problem, we introduce another paraproduct n«(-)) tailor-made to deal 
with that problem, and intertwined to nt,(-) via ^ that is 

on^ = 

so n is formally the 11 operator seen in a different basis 

n = ^"^ono^. 

We show in Section 14.11 that 11 and IT have the same analytic properties. In partic¬ 
ular, if / £ with —2 < a < 0, the Schauder estimate proved in proposition 

[ini shows that ^/) is an element of the parabolic Holder space In the 

end, we shall be working with an ansatz for the solution space of the 3-dimensional 
(PAM) equation given by distributions/functions of the form 

^ + (•••)■ 

The introduction of semigroup methods for the definition and study of para- 
products is relatively new; we refer the reader to different recent works where such 
paraproducts have been used and studied mEKIlElIl]. 


4.1. Intertwined paraproducts We introduce in this section a pair of intertwined 

paraproducts that will be used to analyze the a priori 
ill-posed terms in the right hand side of the parabolic Anderson model equation 
and multiplicative Burgers system in the next section. We follow here for that 
purpose the semigroup approach developed in [3], based on the pointwise Calderon 
reproducing formula 

/= f QfVy TT’fV, 

Jo ^ 

where / is a bounded and continuous function. This formula says nothing else than 
the fact that 

lim = Id. 

tio 

(This is a direct consequence of the fact that the operator ip* tends to the identity 
operator, since ip has unit integral.) We can thus write formally for two continuous 
and bounded functions /, g 


fg = lim Pf) (pi'V . = - £ tdt (pfV • pf^^) } j + A-i(/,5) 

= r (iPfV • Qf^g) + j + A_i(/,g), 

(4.1) 


where 


A_i(/,ff) ■.= VP 


stands for the “low-frequency part” of the product of / and g. This decomposition 
corresponds to an extension of Bony’s well-known paraproduct decomposition m 
to our setting given by a semigroup. 


The integral exponent b has not been chosen so far. Choose it here even and no 
smaller than 6. Using iteratively the Leibniz rule for the differentiation operators 
Vi or dr, generically denoted D, 

D{(pl)(l)2 = D{(pi ■ <t)2) - 01 • T)(02), 
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we see that f ■ 'Pt^^9^ can be decomposed as a finite sum of terms taking 

the form 

where e StGCz and the tuples I, J and integers k,£ satisfy the constraint 

Denote by Ih the set of all such (I, J, k,£). We then have the identity 

r pf> (efV ■ d%) f - S “fi C -"fiT- 


ib 


for some coefficients %’/■ Similarly, we have 


r 


k,e ■ 
Qib) 


dt 


Y.i’ii [‘slh/.s) 

I. J» 


dt 

T’ 


with of the form 

for some coefficients bjj'j,. So we have at the end the decomposition 

f9 = Yi “m + K/idJ)) j + Ti ^k'iU.g) J, 


which leads us to the following definition. 

Definition. Given / 6 Use(o 1 ) 9 ^ L^{A4), we define the paraproduct ng^^(/) 

by the formula 


<>(/):-["{ 2 4’i4'i{/.9)+ E 

and the resonant term Il^’^\f,g) by the formula 
'1 


dt 


y, 


° T,-hi 4 -i.=ih-!-<'=* J 


dt 

t 


ib-,^+k=^+e=i 


With these notations, Calderon’s formula becomes 

fg = nf (/) + uf\g) + n(')(/,5) + A_i(/, 5) 
with the “low-frequency part” 

A_i(/,5) :=V? 

If b is chosen large enough, then all of the operators involved in paraproducts and 
resonant term have a kernel pointwisely bounded by a kernel Qt at the right scaling. 
Moreover, 
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(a) the paraproduct term ng^^(/) is a finite linear combination of operators of 
the form 



with Q\Q‘^e StGCi, and e StGCP’^^]. 

(b) the resonant term (/, g) is a finite linear combination of operators of the 
form 



with Q},Q^e StGCi and e StGC^’^^]. 

Note that since the operators Q* and V} are of the type Q!f \ or a they 

can easily be composed on the left with another operator this will simplify 

the analysis of the paraproduct and resonant terms in the parabolic Holder spaces. 
Note also that 1) = 0, and that we have the identity 

<^(/) = /-pfVf/, 


as a consequence of our choice of the normalizing constant. Therefore the paraprod¬ 
uct with the constant function 1 is equal to the identity operator, up to the strongly 
regularizing operator 


One can prove the following continuity estimates in exactly the same way as in 
[3]. Note first that if u}i,u }2 are two space-time weights, then ui := uiU 2 is also a 
space-time weight. 


Proposition 13. Let uji,u !2 be two space-time weights, and set u := ojiuj 2 - 

(a) For every q;,/ 3 e M and every positive regularity exponent 7 , we have 

\\^-i{f,9)\\c-f ^ ll/llc“ ll5llc/5 

for every / 6 and 5 6 C^^- 

(b) For every a e (—3, 3) and / G , we have 


n(')(/) 


< 


^2^9\\Jf\\cs^ 


C a 

f,! 


for every g e L^{ui 2 ), and 

n(')(/) 


■>a+l3 


< 


bWrP ll/ilcs, 


^^2 


for every g e w'th /3<Oanda-|-/3G (—3, 3). 


(c) For every a, (5 e (— 00 ,3) with a -F /3 > 0, we have the continuity estimate 


n(')(/,5) 


1+/3 


< 


I/lies, II5IIc/5 




for every / G and 9 ^ ■ 

The range (—3,3) for a (or a-F/3) is due to the fact that all the operators involving 
a cancellation used in this estimate satisfy a cancellation of order at least i^-F 10 > 3. 
We simply write 3 in the above statement, which will be sufficient for our purpose. 
We proved similar regularity estimates for the paraproduct introduced in [3], with 
a range for a limited to (—2,1). This difference reflects the fact that the class 
of operators L considered in [3], characterized by the first order carre du champ 
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operators, is more general than the class of Hormander form operators considered 
in the present work, and allows only for a first order calculus. 


These regularity estimates can be refined if one uses the specific weights w and 
PaW introduced in Subsection 13.41 


Proposition 14. For every a e (—3, 3) and a,e e (0,1) with a — a — e e (—3, 3) and 
/ e Cp ^, we have 
• for every 


nW(/) 


ffLII/lb ; 


for every g e with /3<OandQ: + /3 — aG (—3, 3) 


n(')(/) 


-fCx.-\-13 —2(ad-e) 




The proof of this result is done along exactly the same lines as the proof of 
Proposition 1131 using as an additional ingredient the elementary Lemma [TTl 


We shall use the above paraproduct in our study of the parabolic Anderson model 
equation, and multiplicative Burgers system, to give sense to the a priori undefined 
products uC, and of a function rt on A4 with a distribution on 

A4, while 2 q; — 2 ^ 0. Our higher order paracontrolled setting is developed for that 
purpose. As said above, and roughly speaking, we shall solve the Anderson equation 

{dr + L)u = uC 

by finding a hxed point to the map $(n) = e~'^uo + ^~^{uC). We would like to 
set for that purpose a setting where the product uC, can be decomposed as a sum of 
the form 

uC = 2nW(y,)+ (•••), 

i=l 

for some remainder term (•••). We would then have 

j^-H<) = 2^-'(nW(yo)+ (•••), 

which we would like to write in the form 

commuting the resolution operator with the paraproduct. The commutation is 
not perfect though and only holds up to a correction term involving the regularizing 
commutator operator ng(-)], whose regularizing effect happens to be too 

limited for our purposes. This motivates us to introduce the following operator. 

Definition. We define a modified paraproduct setting 

Uf\f) (T/)). 

The next proposition shows that if one chooses the parameters that appears 
in the reference kernels Qt, and the exponent h that appears in the definition of the 
paraproduct, both large enough, then the modified paraproduct has the same 

algebraic/analytic properties as 










31 


Proposition 15. If the ambient space M is bounded, then for a large enough choice 
of constants and 6, the modified paraproduct ^g{f) is a finite linear combination of 
operators of the form 



with e Gct-^ ^ ^ 

If the space M is unbounded, then the result still holds on the parabolic space [0, T] x 
M for every T > 0, with implicit constants depending on T. 

The operators Qj that appears in the decomposition of ng(f) are elements of 
StGC^*^’^^], while the operators Qj that appear in the decomposition of ng(/) are 
mere elements of GCs'^. 

Proof — Given the structure of as a sum of terms of the form 

with e StGC and Q^, e StGC4, it suffices to look at 

We have e StGC^'^’^^], and it is easy to check that also belongs 

to StGC4+^ StGC4. Insofar as 


we are left with proving that the family Q} := belongs to 

GCs”^, with Q} essentially given here by 

Q\ = 

with -^ + /c > |. Note in particular that we have either |/| ^ | or A; ^ We 
check in the first two steps of the proof that Q 6 G in both cases provided h is 
chosen big enough. The third step is dedicated to proving that e GC4~ . 


Step 1. Assume here that |/| ^ The kernel AT of o (t is given by 

coo 

A:((x,T),(y,fT)) = J - A) (4.2) 

So by the Gaussian estimates of the operator tV at scale max(t, A— 

and since |/| ^ we deduce that 


PI {x,y) 


< 


< 


< 


t 


t “t" A — (J 

t 

t “t" A — (T 

t 

t “t" A — (T 


gt+\-aix,y) 


b _1/ 

8 2 


n{B{x, Vi)) M 1 + 


d{x,yf 


2 \ -P 


l-v-b 


n{B{x,Vi)) ^ ( 1 + 


t X — (J ^ 
d{x,y)'^^ 


t 
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ifb is chosen large enough for | —to be non-negative. Using the smoothness 


of y? we then deduce that 


K{{x,T),{y,a)) 


is bounded above by 


n[B[x, Vt)) M 1 + 


d{x,yf 


2 \ 


f 


< 


t ia{B{x,\^)) 


d{x,yf 

t 


2 \ -B 


t “t" A — (7 


r — a\ 


l-v-U 




r — AA dX 

~~r j 


So we get the upper bound 
-fA((x,r),(y,cj)) <v{BM{{x,T),y/t)) 


-1 


d{x,yf 


\T — a\ 


-t-i 


(4.3) 


If d{x,y) ^ 1, this is exactly the desired estimate. If d{x,y) ^ 1 and one works 
on a hnite time interval [0,T] then we keep the information that |A — (t| ^ T 
and so the exponentially decreasing term in the Gaussian kernel on the spatial 
variable allows us to keep in all the previous computations an extra coefficient 
of the form 

- d(x,y)^ 

y{BM{x, i+r 

which is exactly the decay required in the dehnition of the class G. 

Step 2. Assume now that /c ^ |. We work with the above formula for the 
kernel K and use the cancellation effect in the time variable by integrating by 
parts in A for transporting the cancellation from time to space variable. So 
starting from formula (j4.2p . the “boundary term” in the integration by parts 


Vi(r - A) 

is vanishing for A ^ oo, and equal to 

<i){x,y){tdrf-^M'r - <y) 

for A = cr. The latter term satisfies estimate (14.3p . So up to a term denoted by 
(/), bounded as desired, we see that K{^{x, t), (y, a)) is equal to 




(/)+r 

Jo 

where we used that by analyticity of L in L^{M) 




Doing k integration by parts provides an identity of the form 

rOO 

d^((x,T),(y,a)) = (/) + 

Jo 


K \i\ , (x,v)ipt(x — X) 


where (>/) stands for a term with (j4.3p as an upper bound. This procedure 
leaves us with a kernel which has an order of cancellation at least | in space; 
we can then repeat the analysis of Step 1 to conclude. 


Step 3. The proof that actually belongs to GCs”^ is very similar, with 
details largely left to the reader. The above two steps make it clear that the 
study of reduces to the study of operators with a form similar to that of the 
elements of StGC^^’^^^. We have provided all the details in Proposition [3] of how 
one can estimate the composition between such operators and obtain an extra 
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factor encoding the cancellation property. The cancellation result on Q} comes 
by combining the arguments of Proposition [3] with the two last steps. 

Let us give some details for the particular case where the family Q belongs to 
StGC“ for some a ^ | — 1 and commutes with this covers in particular 

the case where Q is built in space only with the operator L with no extra Vi 
involved. Let us then take s,t e (0,1) and consider the kernel of the operator 
QtQl- Note hrst that 

qIq: = {qIq:) o 
1 - 1 

Since Q e we know that QfQ* is an operator with a kernel with decay at 

b 

scale (t + s )2 with an extra factor ^ ^ We may also consider that 

±_ 

2'^: = ((rriyi)" ^ 


for some operator having |-order of cancellation and a kernel with decay 
at scale ^/s~^Vt. So by what we did in the two hrst steps we also obtain that 
Qt+si^ + has a kernel with decay at scale {t + s)^, for a large enough 

choice of b. (Indeed, note that is very similar to the operators studied in the 
two hrst steps: easily analyzed as a function of the space-variable, while, as far 
as the time-variable is concerned, the composition of convolution preserves the 
main properties needed on the functions ~ vanishing moments.) At the end, we 
conclude that 



iit + sr) 


— —1 
16 ^ 


Q 


2 

£+s 


with Qf_|_s having fast decreasing kernel at scale (s -I-1) 2 . That concludes the 
fact that Q} e GCs"^. 


[> 


The following continuity estimate is then a direct consequence of Proposition [T21 
since the latter implies that we can reproduce the same argument as for the standard 
paraproduct in Proposition [TH 

Proposition 16. For every a e (—3, 3) and a, e G (0,1) with a — a — e e (—3, 3) and 
/ e Cp^, we have 

||nf(/) 

for every g e L“. 

Last, note the normalization identity 

ni(/) = / (Jf/) 

for every distribution f e it reduces to 

ni(/) = /-7’fVf^(/) 

if /| ^ = 0. (Use here the support condition on (p in the definition of V.) Let us 

also point out here the strongly regularizing effect of the two operators and 


< 


||ro ^g\\ 


ll/ll 
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^ denoted by A below, that satisfy the continuity estimate 



for any a, /3 e (—3, 3) and any space-time weight u. 

We shall fix from now on the parameters b and ii, large enough for the above 
result to hold true. 

Remark 17. The previous Proposition is very interesting because of the following ob¬ 
servation: the time-space paraproducts IT are defined in terms of parabolic cancellations 
and so do not differentiate the space and the time. Consequently, it is not clear if the 
time-space paraproducts IT may be bounded on for some a < 0 (with or without 

weights). Such property would be very useful since the paracontrolled calculus (as shown 
later in the study of (PAM) for instance) needs to estimate the composition of (the 
resolution of heat equation) with the paraproduct. However, following the definition of 
the paraproduct we have for / e and g ^ 

(/) = 

So if / e for some a e (—2,0) then Schauder estimates imply that f e 

and we may then use the boundedness on Holder spaces of the modified paraproduct 

In conclusion, these new space-time paraproducts seem to be very natural for the 
paracontrolled calculus. They allow us to get around a commutation between the initial 
paraproduct and the resolution operator (which could be a limitation for a higher 
order paracontrolled calculus) and fits exactly in what paracontrolled calculus requires 
to solve singular PDEs, modelled on the heat equation. 


4.2. Commutators and correctors We state and prove in this section two continuity 

estimates that will be useful in our study of the 
3-dimensional parabolic Anderson model equation and Burgers system in Section [5l 

Definition 18. Let us introduce the following a priori unbounded trilinear operators on 
5'. Set 

R(f,g,u) :=nW(nf(/))-ng(/), 

and define the corrector 

C(f,g,u) :=nW (nf (/),u) -<7n(')(/,u). 

This corrector was introduced by Gubinelli, Imkeller and Perkowski in |22] under 
the name of commutator. We prove in the remainder of this section that these 
operators have good continuity properties in some weighted parabolic Holder spaces. 

Proposition 19. Given some space-time weights wi,W 2 ,W 3 , set co := Let 

a,^,^ be Holder regularity exponents with a e (—3,3), (3 e (0,1) and 7 e (-3,0]. 
Then if 6 := a + /3 + j e (—3, 3) with a + /3 < 3, we have 

||^(/,5,^^)|lc^ ^ ll/llcs, bWc^^ Mcz^, (4.4) 

for every / e g e u e CZ^', so the modified commutator defines a trilinear 

continuous map from x 0^2 x to C£. 
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Proof — Recall that is given by a finite sum of operators of the form 

) :-£ci'(Q?(-)n'(9)) J, 

where belong at least to StGC^. We describe similarly as a hnite 

sum of operators of the form 

4 (-)at 

Thus, we need to study a generic modified commutator 

and introduce for that purpose the intermediate quantity 

£{f,9,u) := Qf{Qt{f)-Vl{g)-VUu)) 

Note here that due to the normalization IIi ~ Id, up to some strongly regular¬ 
izing operator, there is no loss of generality in assuming that 

Ql-Q^t j = £ Ql' Qtj = Id. (4.5) 


Step 1. Study of [Ag{f)) — £{f,g,u). We shall use a family Q in StGC“, 
for some a > |(I|, to control the Holder norm of that quantity. By definition, 
and using the normalization (14.5p . the quantity Qr(^Au {Ag{f)) —£{f,g,u)^ is, 
for every re (0,1), equal to 


£ £ Q.Qr{Q"Qt*( 2 ?(/)n'(ff)) ■ -pU^)} ^ 
= ££ QrQl‘{QtQl'{QUf){'Pti9)-'Pli9))) ■ 



QrQl-{Qt{f)-Vl{g)-V^{u)) 



dsdt 
st ’ 


ds 

s 


where in the last line the variable of Vgid) Qs*^ frozen 

through the action of Then using that g e with (3 e (0,1), we know 

by Proposition [3 that we have, for r ^ u. 


cv 2 (x,t) ^ (Pjg)(x,T) - (Plg)(y,a) < (s + t + p((x,t), (y,cr)) 


2 gC(i(a;,i/)j|^|| 


Note that it follows from equation (13.41) that the kernel of QgQj* is pointwise 
bounded by Qt+s, and allowing different constants in the definition of the Gauss¬ 
ian kernel Q, we have 


St+s((x,T),(y,a)) (s + t + d(x,y)^)^ < (g + t)2 gt+g((x,T), (y, a)). (4.6) 

So using Lemma [Hand the cancellation property of the operators Q at an order 
no less than a (resp. 3) for Q (resp. the other collections Q*), we deduce that 


w ^Qr(^A^(Al^(f))-£(f,g,u)^ 


00 


< 


l/llc“ WaWrpMcz 



iu+tt) 


2 ., . I 1 dsdt 

t 2 {s + t) 2 s 2 


st 
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where we used that 7 is negative to control Vg{u). The integral over t e (0,1) 
can be computed since a > —3 and a + /? < 3, and we have 

u;-^Qr(^A^Al(f)) - £(f,g,u)^ 



uniformly in r e (0,1) because |o| > S. That concludes the estimate for the 
high frequency part. We repeat the same reasoning for the low-frequency part 
by replacing Qr with Qi and conclude that 


Al {Alif)) -£{f,g,u) 


< 


Ct 


llr^ • 


Step 2. Study of A^g — £{f,g,u). This term is simpler than that of Step 1 
and can be treated similarly. Note that Qr(^Ag {Au{f)) — £{f,g,u)^ is equal, 
for every r e ( 0 , 1 ), to 

£ QrQ!-(Qt(f)r^(ug)) ^ - £ QrQ!-(Qt(f) ■ Vl{g) ■ Vl{u)) ^ 

= £ QrQl'(Qt{f){Vl{ug) - Vl{g) ■ Vl{u))) 

Now note that since g e with /3 e (0,1), we know by Proposition 0 for 

T ^ a, 


W 2 (x, r) ^ g{x, t) - {Vlg) (y, a) 

< a; 2 (x, t)“^ y(x, r) - y(y, a) + C 02 {x, t)~^ g{y, a) - {Vjg) {y, a) 

,2\ 2 


< 


(^s + t + p{{x, t), (y, a)) \\g\\^p ^. 


Then the same proof as in Step 1 can be repeated. 


[> 


As far as the continuity properties of the corrector 

c(/,y,u) = n(')(nf (/),«) -ynW(/,n) 

are concerned, the next result was proved in an unweighted setting in [U Proposition 
3.6] for a space version of the paraproduct 11; elementary changes in the proof give 
the following space-time weighted counterpart. 

Proposition 20. Given space-time weights 001 , 002 , 003 , set 00 := 001002003 . Let a,(3,j 
be Holder regularity exponents with a e (—3,3),/? e (0,1) and 7 e (— 00 ,3]. Set 
5 := (a -I- /3) A 3 -I- 7 . If 

0 <a-|-/ 3 -f 7 <l and a -I- 7 < 0 

then the corrector C is a continuous trilinear map from x x CZ^ to C^. 
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5 _ 

Anderson and Burgers equations in a 3-dimensional background 

We are now ready to start our study of the parabolic Anderson model equation 

{dt + L)u = uC, 

and the multiplicative Burgers system 

{dt + L)u + {u ■ V)u = 

in a 3-dimensional manifold, using the above tools. Here for Burgers system, we 
consider a collection of three operators V := (Vi,V 2 jV 3 ), so^o = 3 here. We shall 
study the (PAM) equation in a possibly unbounded manifold, using weighted Holder 
spaces, while we shall be working in a bounded setting for the Burgers equation, as 
its quadratic term does not preserve any ’obvious’ weighted space. 


5.1. Getting solutions for the (PAM) equation Let us take the freedom to assume for 

the moment that the noise ( in the 
above equations is not necessarily as irregular as white noise. We shall hx from now 
on a hnite positive time horizon T. Recall the elementary result on paracontrolled 
distributions u with derivative u stated in section[2l such distributions are of the form 
u = e~^vi, for some more regular factor vi. This is indeed what happens formally for 
any solution to the (PAM) equation, since uC, = nu(((), up to some smoother term, 

and (^u{C)^ = ^u{-^~^C), up to some more regular remainder. Elaborating 
formally on this remark leads to the introduction of the following distributions, and 
the choice of representation for a solution of the (PAM) equation adopted below in 
Proposition [21] 

For a continuous function ( in , and 1 ^ i ^ 3, dehne recursively the following 
reference distributions/functions 


with 

£o £o 

Fi:=C, H 3 := 2 ^*(^i)^i(^2), (5.1) 

and define 

to to 

(*) := -2Y,y{Zi)Vi{Z^), W 2 :=-^Vi(Z2f 


as well as for j 6 {1, ..,.^ 0 } 


Wi : = 


^n(^^(^Vi(Zi),v,^-\VjZi)y 


Indeed in the term (*), only the resonant parts in the products have to be dehned, 
since the parapoducts always make sense, so we focus on the resonant part of (*) 


to 

Wi ■.= -2'^n(^)(v,(Z3),Vi(Zi)). 
i=l 

Defining 

• the y/s as elements of c 
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• the distributions as element of L^Cp^ 

• the quantities W 2 as elements of L^Cp^~^, 

for some 1/3 < a < 1/2 and a > 0, when (/ is a space white noise, is the object of 
the renormalisation step, which shall be done elsewhere. These conditions ensure, 
by Schauder estimates, Proposition [9l that Zi is in the parabolic Holder space C*". 
Note that assuming Wi is an element of L'^Cp^~^ ensures that (★) is an element of 
L^Cp~^. There is a clear correspondence between the terms defined here and those 
appearing in Hairer and subsequent analyses of the KPZ equation; see e.g. 1201 [211 
I28| . In a simplified setting where the vector fields Vj are constant and correspond 
to the derivation operator in the direction of the vector of the canonical basis, 
the above terms correspond to 

IT/ = {dZi) • d^-\djZi), Z 2 = ^-\idz,f), 1 T 2 = {dz^f 

Zs = ^-\{dZi){dZ2)), ITi = {dZi){dZs). 


Set 


Z := Zi + Z2 + Z^ =: Zi + Z. 


Proposition 21. The function u is a formal solution of the (PAM) equation if and only 
if the function 

V := e~^u 


is a solution of the equation 


4 

^v = -Uv + 2Y,V^{Z)Vi{v), (5.2) 

with the same initial condition as u at time 0. The letter U stands here for ITi + IT 2 + IT 3 
for an explicit distribution IT 3 in 

We explicitly single out IT 2 here, and not IT 3 , even though they both belong to 
the same space, because IT 2 will later have to be renormalized while IT 3 will be 
well-defined as soon as the other quantities Yi, ITi,... will be well-defined. 

Proof - Observe that 


drU = ydrV + vdrZi + vdrZ 

and using the Leibniz rule on I/’s 
fo 


Lu = 2 y^{Zfv - Vl{Z)v - 2T(Z)T(u) - 

= ^vLZ + Lv- 2 Vi{Zfv - 2T(Z)T(u)^ 

= ivLZi + vLZ + Lu - u 2 Vi[Zf - 2T(Z)T(u) | . 


2=1 
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Due to the definition of l^’s, we have some telescoping property: 
^Z-Y,V^[Zf 

= + Is) - 2 [V^iZl) + V,iZ2) + ViiZs)) 

i=l 

4 3 

= h 2 + Ds - ES V^{Zj)Vi{Zk) 

i=lj,k=l 
io 3 

= Vhi + Vh 2 - E E V,{Zj)V^{Zk). 

j,k^2 
j + k'^5 

Since we assume that Zj e it follows that Vi{Zj) e and Vi{Zk) 6 

L’^Cp^~^. Given that j + k ^ 5 and a e (1/3,1/2), at least one of the two 
numbers {ja — 1 ) and {ka — 1) is positive and the other not smaller than 2 q: — 1 . 
So 


^0 

U :=^Z -Y, £Wi + W2 + (5.3) 


and the result follows. 


[> 


We solve (15.2p using paracontrolled calculus instead of solving directly (PAM). 
Definition 22. Given | < /3 < a < ^ and a time-independent distribution ( e ^ 

(PAM)-enhancement of C is a tuple Q : = (C,l2,l3,Wi,W2,(W2^),), with 


Yk e 


and 


Wi,W2,Wi£ 


cc/~<2a—l 


So the space of (PAM)-enhanced distributions C, for the (PAM) equation is 
here simply the product space 


00/o2a— ®(^ 0 + 2 ) 


qr" X n ircrX 




5 . 1 . 1 . The paracontrolled approach The study of singular PDEs, such as the An¬ 
derson and Burgers equations or (15.2p , is a four 
step process from a paracontrolled point of view. Let us sketch it for equation (|5.2I) 
as an example. 

(a) Set yourself an ansatz for the solution space, in the form of a 
Banach space of paracontrolled distributions/functions. 


Given ^</3<a<^,we choose here to work with functions v paracontrolled by 
the collection (Vi(Zi)) | , that is with v of the form 

^0 

v = (5.4) 

i=l 

for a remainder v'^ e ^ and Vi e C^. We refer the reader to Subsection 13.41 for 

the introduction of weights pa and w. Note that we use the IT paraproduct and not 
the n paraproduct. We turn the solution space 

‘5a,/ 3 (^ := |(u; ui,... , u**) satisfying the above relationsj 

into a Banach space by defining its norm as 

4 

||('(;;ui,... := ||'(;‘'||^ofa+/3 + ^ \\vi\\(,g- (5-5) 

2 = 1 

(b) Recast the equation as a fixed point problem for a map <f> from the 
solution space to itself. 

This is where we use the continuity properties of the corrector and different para- 
products. In the specific situations of equation (15.2h . given {v,vi,... ,V£q‘,v^) in the 
solution space Sa ,/3 , one sets 

y = ^(^ - Uv + 2Y,V^{Z)V^{v)j 
2^1 

and shows that it has a decomposition (y; yi,..., yt*) of the form (15.4p . This 
is where we need all the extra information contained in Then, given an initial 
data Vo e , the application 7 : (r,x) e~'^^{vo){x), belongs to and 

satisfies 

= 0, 7r=o = ^^0- 

We define a continuous map 4* from the solution space to itself setting 

4> := ^ (y-h 7 ; yi,..., y^^; y** -^ 7 ). 

(c) Prove that $ is a contraction of the solution space. 

Recall a parameter k > 1 appears in the definition of the special weight w. We 
shall see below that the function y** satisfies the estimate 

j|yt*j|^nf «+/3 ^ K~^\\{v;vi,... 

for some e > 0, and that (yi,..., y^g) depends only on v and not on ui,..., and 
. These facts provide a quick proof that 4> o is a contraction of the solution space 
‘5a,y • Indeed, given (u; ui,..., vg ^; ut*) in 5 a ,/3 , set 

{z -t- 7 ; 2 ;i,.. ■,zg^-,z^ + 7 ) := ^°‘^{v,vi ,... e ‘ 5 a,/ 3 (^. 

We know that 
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The paracontrolled structure (15.41) of y and Schauder estimates also give 

^0 

WvWr^+P ^ i^~^\\yi\\c^+°‘ 

'"-rop-a ^ 

i=i 

So we conclude that y is controlled with a small bound. Since (^i,.., zi^) will be given 
by y, we then obtain that {z;zi,..., zi^; z^) will be controlled in Sa^/ 3 [C ) with small 
norms, relatively to the initial {v;vi,... ,ViQ;v^), so $ o <f> will indeed be a small 
perturbation of the map {v;vi,..., V£q;v'^) ( 7 ; 0,... , 0; 7 ). Then it is standard 

that if K~^ is small enough, that is k is large enough, then we can apply the hxed 
point theorem to and conclude that it has a unique fixed point in the solution 
space Sa,j 3 (C); the same concolusion for < 1 > follows as a consequence. 

(d) Renormalisation step. 

The dehning relations for Zj found in step (b) actually involve some terms that 
cannot be dehned by purely analytical means when ^ is a white noise, but which 
make perfect sense for a regularized version of C- Their proper dehnition requires 
a renormalisation procedure that consists in dehning them as limits in probability, 
in some parabolic Holder spaces, of suitably modihed versions of their regularized 
versions (with in place of ("), which essentially amounts in the present setting to 
adding to them some deterministic functions or constants. (This may be trickier 
in other situations as the theory of regularity structures makes it clear.) Given the 
inductive construction of the Zj, this renormalisation step also needs to be done 
inductively. At e hxed, this addition of deterministic quantities in the dehning 
relations for Zj dehnes another map from the solution space to itself that can 
eventually be equivalent to consider a renormalised equation with noise with 
e-dependent terms added in the equation, when compared to the initial equation. 
Write for its solution. In the end, we get, from the continuity of hxed points of 
parameter-dependent uniformly contracting maps, a statement of the form: Let 
stand for the map constructed by taking as reference distributions/functions Zi the 
limits, in probability, of their renormalised versions. Then the functions converge 
in probability to the solution u of the fixed point problem of the map $. 

We shall do here the hrst three steps of the analysis for both the Anderson and 
Burgers equations, leaving the probabilistic work needed to complete the renormal¬ 
isation step to another work; we shall nonetheless give in Section [ 6 ] some hints as to 
what is going on. 


5.2. The deterministic PAM equation Given what was said in the preceding section, 

the main work for solving the (PAM) equa¬ 
tion consists in proving the following result. 

Theorem 23. Let ^ < a < ^ be given. Choose j 5 < a, the positive parameter a in the 
weight Pa, and e > 0, such that 


2a + > 1 and 8 (a -I- e) ^ a — /3. 
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Given an enhanced distribution C, one can extend the product operation 

^0 

V e C®(M) -Uv + 2 2 Vi{Z)Vi{v) 

i=l 

to the space Sa,/3{C) into an operation v i—> —Uv + 22^2=1 Vi[Z)Vi{v), so that setting 


y := 


l-o 

Uv + 2Y,Vi{Z)Vi{v)\, 

i—1 


and yi := 2vVi{Z^) + 2Vi{v), there exists 6 such (r/;r/i, 

element of the solution space Sa,i3{C), and 


{y;yi,---,y£o\y^) 


< 


Q!,/3 


\yH i+c+p < K 


{v,vi,.. .,Vig;v^) 


a,0 

O', 13 


yio;y^) is an 


(5.6) 


Proof — First, we note that since v satisfies the ansatz (15.41) and 2a < a — /3, we 
know from Schauder estimates that 


V e C 


l + Q 
-^Pa 


n 


'■'■ujp-a ■ 


Step 1. We first consider the part Uv where we recall that U = Wi + W 2 + VF 3 
for some VF 3 e L^Cp'^~^. Using the paraproduct algorithm, one gets 


W 37 ; = u^^^iv) + n(')(W3) + u^^\v,W3). 

By the boundedness of paraproducts. Proposition [131 and Schauder estimates. 
Proposition [T2l we get 


with 


nU^ 3 (^;) 6 so {v) e clfl cz ci+p^_+f 




^l + a+13 


n ^\\v 


\^i+p 


^ — a 

since 2e + 2a < a — /3 and a < 1. For the 
wifVi Prnr\n«itinn II .‘^1 vipIHcj 


rpQriTiPTif ■nprf 


U^^\v,Ws)eCl^+^ so 




with 




^l+a.+13 


< K 


V\\^l+I 3 . 


For the second paraproduct, we use the modified paraproduct and its bound¬ 
edness, Proposition [161 to have .if'“^ni^^(lU 3 ) = hence since 

J^-ilU 3 e Ci+ 2 " we have ^-^U^v\W 3 ) e Ci^p‘lt^ with 




^1-f q; + ,S 




"U It . 

"^zupa 


since 4(a -I- e) ^ a — /3. So we have ^ ^(Wsu) e , with an acceptable 

bound. 

The term IF 2 is an element of L^Cp^~^, so using the same reasoning yields that 
J^“^(VF 2 u) e with an acceptable bound. 
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' oo/^a —1 

'-'Pa 

two previous terms. Recall its definition 


The term Wi is an element of so it is really more singular than the 


Wi = -2£Vi(Zi)V,(Z3) 


i=l 


with Vi{Z 3 ) in since Z 3 is an element of So Wi is in Cp^ and since 

veClt^, we have 

^Pa ‘ 


and nW(lTi,u)eC: 


Using Schauder estimates one obtains 


2a 

'^P2a 


Ql + a + fS 


(^l + a + p 


< K ^llulTl + a. 


It remains us to study the paraproduct term 

nW(IUi) = 1 + 11 + 


with 


I:. -2|;n('')(<>2^,+.(Zi))) 

i=l 

eo 

II := -2 2n(')(n(')(u(Z3),u(Zi)) 

^0 

Ml:. -22n(‘)(niJ>^_,(v;(Z3))). 

i —1 

By easy considerations on paraproducts, the third term III belongs to 
and ^~^{\\\) e with acceptable bounds, because Z 3 is an element of 

Cp+". Moreover, since we assume that VUi = (Uj(Z 3 ), Vi(-2’i)) is an 

element of L®Cp"“^, the second term II also satisfies ^~^{\\) 6 Crop'^^^. Using 
the regularity of u 6 Crop^ ^mpa aad Proposition [19] for the commutation 
property, we dednce that 


74q!-2 

-'rop3a 


i=l 

and consequently 

£•-‘(1) € -2 2 n<« ,, [^-V,(Z 0 ] + 


2 = 1 


with an acceptable bound for the remainder since 8 (a + e) + 1 < 3a — /3. 

At the end, we have obtained that 

Ji-‘(uv) 6 {22 + ci+“y}. 

2^1 

which proves that ^~^(Uv) is paracontrolled by the collection j 

and the remainder has a bound controlled by k~^. 

Step 2. Let now focns on the term Xlili Fix an index i and write 

U(z)u(u) = n^^)^)(U(z)) +n(')(U(z),u(u)). 
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The second term is of regularity 2a — 1 and using the modified paraproduct, 
Schauder estimate and the fact that we have v e Cii", we see that 

zupa ' 


^-1 


n 


( 6 ) 


(^*(^)) 


= n 


( 6 ) 

Vi(^) 




We proceed as follows to study the resonant part. First, since a > 1/3, we have 


^0 

1=1 


Consider the modified resonant part 

Uf\f,g) :=n(')(/,F.3) 


and the corresponding corrector 

Ciif,g,h) :=TTf (nf (/),h) (/,h). 

Then since in the study of the resonant part and the commutator, we can change 
the localization operators, so we can integrate an extra Vi operator, we get 
boundedness of from C" x to as soon as a + /3 — 1 > 0, and 

boundedness of the corrector Cj from C“ x X into C'^+^+'y 1 as soon as 
a + /3 + 7 — 1>0, proceeding exactly in the same way as above for and C. 
Using this commutator, we see that Xji=i is an element of the 

space 

£o ^0 £o 

Y,Y^v,.u(^\Vi{Z,),V^-\V,Z^)) + Y,a{V{Z,),v,,Vi^-\V,Z,)) 

i=lj=l i=l 

that is an element of 


1=1 

since W 2 e and 2a + /3 > 1. In the end, we conclude that 


y.= ^-^ -Uv + 2'^ViZ)Viiv) 


as expected. Observe that 1 /(^ 3 ) is of parabolic regularity (3a 
and Vi{v) belong to . 


1 ), so vVi{Zs) 
> 


We can then apply the contraction principle, such as explaned above in Step (c) 
in section 15.1.11 

Given uq £ , write ) foi those tuples (u; ui,..., V£q', v'^) in (C) 

with U|i-=o = vq. As the function 7 := (x,r) 1-^ (uo)(x) belongs to 

and is the solution of the equation 

(d.r + T)( 7 ) = 0 , 7^=0 = 110, 

we define a map from (C) to itself setting 

$(u;ui,.. .,vy,v^) = {y + 7;2uUi(Z3) + 2Ui(u),... ,2vVig{Z3) + 2Veg{v)]y^ + 7), 
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with 

y:=^-^[-Uv + 2 Y, Vi{Z)Vi{v)), 
i—1 

and yt* given by the previous theorem. Note that the map ^ depends continuously 
on the enhanced distribution the next global in time well-posedness result is then 
a direct consequence of Theorem 1231 

Theorem 24. Let us work under assumption (A), and let | < a < ^ be given. Choose 
/3 < a, the positive parameter a in the weight pa, and e > 0, such that 

2a + j 3 > 1 and 8 (a + e) ^ a — /3. 

Then, one can choose a positive parameter k, in the definition of the special weight zu, 
large enough to have the following conclusion. Given vq e the map has a 

unique fixed point (u,ui, ..,V£g,v^) in it depends continuously on the enhanced 

distribution (, and satisfies the identity Vi = 2vVi(Z3) + 2Vi{v) for i = 1 ,..,£q. This 
distribution is the solution of the singular PDE 

^0 

^v = -Uv + 2Y, Vi{Z)Vi{v) ( 5 . 7 ) 

with U|,-=o = Vq. The function u = e^v is then the unique solution of the (PAM) 
equation with initial data uq. 

If the ambient space M is bounded, then we do not have to take care of the infinity 
in the space variable, and one can prove a global (in time) result by considering the 
weight ■a 7 (x, r) = e'^'^ with a large enough parameter k. 

5 . 3 . The stochastic PAM equation Recall the time-independent white noise over the 

measure space (M, p) is the centered Gaussian 
process ^ indexed by T^(/r), with covariance 

= J fix)p{dx). 

It can be proved [3] to have a modification with values in the spatial Holder space 

Cp^ , for all positive constants e and a, where u is the Ahlfors dimension of 
(M, d, p) - its dimension in our Riemannian setting. We take z/ = 3 here. We still 
denote this modification by the same letter The study of the stochastic singular 
PDE of Anderson 

can be done in the present setting. This requires a renormalisation step needed 
to show that the quantities E = can be dehned as elements of suitable 

functional spaces, as limits in probability of distributions of the form — A^, where 
S'" is given by formula (|5.1I) with C = the regularized version of the noise 

via the semigroup, and are some deterministic functions. This renormalisation 
step is not done here; Section [U] gives however a flavour of what is involved in 
this process in the present setting. Note that the two dimensional setting was 
studied in depth in [3], with spatial paraproducts used there instead of space-time 
paraproducts. We then formulate this renormalisation step as an assumption in the 
present work. Recall the definition of Z 2 , Wi , W 2 , W 2 , P 3 and Z 3 given in Section 

EH 
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Assumption (B) - Renormalisation. Let ^ stand for white noise on M , and for e > 0, 
denote by := its regularized version, and by the distributions corresponding 

to H = Yi,Zi,y 2)-^2 that one obtains by replacing C, by 

(a) There exists a family (Af)o<£<=:i of deterministic functions such that Y 2 — Af is 
e-uniformly bounded and converging in CrCp^ for every a e (0,1) and any 
a < 1 / 2 . 

(b) Use the upper e-exponent in H = Z 2 ,Wi,W 2 ,W 2 ,Ys, to denote the quan¬ 
tities that one obtains by replacing Zi by Zf, and Y 2 by Y 2 — Af. For any 
a < 1 / 2 , 

• the distributions Z|, y/, are e-uniformly bounded and converging in 

CrCp^"^, respectively and CtCp^~^, for every a e (0,1); 

• there exists deterministic functions Al 1 and Al 2 such that the distributions 
Wl — All and fy| “ -^ 22 ’ are e-uniformly bounded and converging in 
CtCp^~^ and respectively CtCp^~^, for every a e (0,1). 

This assumption about the renormalisation process for the above quantities prac- 
tiaclly means that one can renormalise the most singular quantity Yl by substract- 
ing an e-dependent deterministic function, and that once this has been done, no 
extra renormalisation is needed for the terms Y^ and W^’^. At the same time, the 
quantities Wf and VFI have to be renormalized and this operation can be done by 
subtracting deterministic functions - essentially their expectation. 

Write Z‘" and for the renormalized versions of Z^ and U^. By tracking in the 
proof of Theorem 1231 the changes induced by such a renormalisation of Y^, Wf and 
W 2 , we see that if (v^;vf ,..., satisfies ansatz (15. 4p with ( = and setting 

:= yf’-' (-([/" + Af - A|,1 - AI 2 X + 2 2 V.(^^mvn) , 

i^l 

then the tuple 

(y^;2v^Vi(Zl) + 2Vi(v^),.. .,2v^V,,(Zl) + 2V,,(v^);y^’» + 7 ) 

also satisfies the ansatz. The renormalisation quantity Af — A| 1 — Al 2 comes from 
the definitions of U and U^, since we have after replacement of Yf by y 2 "^ = Y 2 — Af 

l/^=W + Yi-^ V,(ZJ)V,(ZI) 

i=l j,k=l 

= y/ + y|-2 2 v,(z^mz^,)-Af 

i=lj,k=l 

i^l j,k^2 
j + k^5 

_ _ €0 3 

= 2 y,(y|)y,(yO-A! + A|,i + A|,2. 

2^1 j,k^2 
j + k^5 
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Theorem 25. Let us work under assumptions (A) and (B), and let ^ < a < ^ be 
given. Choose f3 < a, the positive parameter o in the weight pa, and e > 0, such that 

2 a + /3 > 1 and 8(a + e) ^ a — /3. 

One can choose a large enough parameter k in the definition of the special weight w 
for the following to hold. There exists a sequence of deterministic functions (-^j)o<£<i 
such that if stands for the solution of the renormalized equation 




to 

- (r + Af - - xy)v^ + 2 2 V,{Tmv^) 

i=l 


v%0) = VO (5.8) 


with initial condition vq e , then converges in probability to the solution 

V e of the Cole-Hopf transformed (PAM) equation (15.7|) constructed from the 

enhancement of the noise given by assumption (B) and Theorem [23l 


By reproducing the calculations of Subsection 15.11 we observe that is solution 
of equation (|5.8h if and only if := is solution of the equation 

= iff — Af + A| + A| 2 )^^) ^*"(0) = Vq. (5-9) 


Theorem 26. Let us work under assumptions (A) and (B), and let ^ < a < ^ be 
given. Choose (3 < a, the positive parameter o in the weight pa, and e > 0, such that 

2 a + /3 > 1 and 8 (a + c) ^ a — /3. 

One can choose a large enough parameter k in the definition of the special weight w 
for the following to hold. There exists a sequence of deterministic functions (-^j)o<£<i 
such that if stands for the solution of the renormalized equation (15.911 with initial 
condition vo e , then converges in probability to the solution u e of the 

(PAM) equation constructed from the enhancement of the noise given by assumption 
(B) and Theorem [231 

This result is coherent with the result of Hairer and Labbe proved in [28]. Indeed, 
in [251 Equation (5.3)], the quantities involving an odd order of noises need no 
renormalisation terms, like us; nor does the term which is of even order but 
involves an extra derivative Vj. We give more insights on the latter term at the end 
of Section and explain why this extra derivative with anti-symmetry properties 
implies that the potential renormalisation term is actually null, as in [28j . 


5.4. The multiplicative Burgers equation We study in this last section the multi¬ 
plicative Burgers system 

{dt + L)u + {u ■ V) u = 

in the same 3-dimensional setting as before with three operators V := (Vi,V 2 ;V 3 ) 
forming an elliptic system. Here the solution u = {vf ,u^) is a function with 
M^-values and {u ■V)u has also 3 coordinates with by definition 

3 

[(u • V) := ^ uWj(tf^). 

To study this equation, we have to make the extra assumption that the ambient 
space M is bounded. Indeed the boundeness of the ambient space is crucial here, as 
using weighted Holder spaces, it would not be clear how to preserve the growth at 
infinity dictated by the weight when dealing with the quadratic nonlinearity. In such 
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a bounded framework, we do not need to use spatial weights and consider instead 
the unweighted Holder spaces - or rather we work for convenience with a weight 
in time 

vj{x,t) := e'^'^. (5.10) 

We stick to the notations of the previous section. The study of Burgers’ system 
requires a larger space of enhanced distributions than the study of the 3-dimensional 
(PAM) equation; the additional components include those quantities that need to 
be renormalised to make sense of the term (tt • V) u, when C, is an element of 
such as space white noise. 

We Hrst rewrite Burgers system in a more convenient form, as we did for the 
(PAM) equation. For each cooordinate exponent j = 1,2,3, we define Za,W^ from 
as above. Then consider a function u : M t—> defined by 

with V : M (—> Then observe that u is formally a solution of 3-dimensional 

Burgers system on M if and only if v is the solution of the system 

3 3 

W = -Wv + 2 2 Vi{Z^)Vi{v^) - 2 {ViV^ + v^ViZ^). (5.11) 

i—1 i—1 

To treat the nonlinearity, we need to introduce another a priori given element in the 
enhancement of the noise Q. Define a 3 x 3 matrix 0 setting formally 

Qij =Yi^^)(^Z{,ViZ{). (5.12) 

Definition 27. Given | < /3 < a < ^ and a time-independent distribution Q e a 

(3d Burgers)-enhancement of c is a tuple C := (c, ^" 2 , ^ 3 , Wi, Ws, (W|),-, o) , with 

Yk e Wi,W 2 ,wi e and 0 e 

So the space of enhanced distributions for the multiplicative Burgers system 
is the product space 

3 

fc =2 

we slightly abuse notations here as the first factors in the above product refer to 
M^-valued distributions/functions, while the last factr has its values in M®. Given 
such an enhanced distribution (, we define the Banach solution space Sa,i3{C ) as in 
Section [5.1.11 replacing the weight pa by the constant 1. Recall the constant k> 1 
appears in the time weight ()5.10p . 

Theorem 28. Let us work under assumption (A), and let ^ < a < ^ be given. Choose 
(3 < a, the positive parameter a in the weight pa, and e > 0, such that 

2a + P > 1 and 6e ^ a — /3. 

Given an enhanced distribution and ve Sa,j 3 {C), the nonlinear term 

3 

[iV(u)]^' := 2 v^e^^{diV^ + v^ViZ^) 

i=l 

is well-defined and there exists some z'^ e with 
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and 




a,/3 


< K 


{v,Vi,..,V3,V^) 




(5.13) 


Proof - We fix a coordinate j = 1,2,3 and have to study 

3 

[Ar(^;)]i = 2 v^e^\ViV^ + v^ViZ^). 
i=l 

The first quantity is sufficiently regular by itself, and we have Zi e C°^, v e 
so for every i = 1,2,3 then 

v^e^^ViV^ e C“ 

licncG 

ZZ-^[v^e^WiV^] e C^+“+^ 
with an acceptable norm (controlled by k~^). 

Let us now focus on the second part v^ViZK Since v e C^°‘, it is very 
regular and the problem only relies on defining the product e^^ViZK We first 
decompose using paraproducts 


e^'ViZ^ = Ii^%{V^Z 


n 


{b) 

ViZi 


{e^') +I[^^\e^\ViZ^). 


The second term is bounded in The last resonant part is studied 

through a paralinearization formula (see [1] and references there for example) 

=n 2 ,(^*) + ( 2 a) 

which implies with a > 1/3 

ylb ;= Ili>’\e^\ViZ^) = (u^^^,{Z^),ViZ^^ + (3a - 1) 

= e^^U^^\z\ViZ^^ + (]3a- 1) 

= e^'U^’’\zi,ViZ() + (3a - 1D = + (3a - 1), 


where we have used the commutator estimates. Since we assume that 0 is 
supposed to be well-defined we conclude to e So we 

observe that 

is well-defined in whose evaluation through A!f~^ is then bounded in 
with acceptable bounds. And since 

this is also controlled in by Schauder estimates and we conclude to 




•71 + 0 ! + ^ 


It remains the quantity with (instead of A*-^). Here we only know that B^^ 
belongs to (and not as for A*-^) but we can take advantage of 

the fact that is a paraproduct. Indeed as before we have 

+ ji^^\b^^ yv^) 

well-controlled in and 
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which is well-controlled in ^ 1 +“+/^ ^ 1^0 to Schauder estimates, Proposition [TOl 
In conclusion, we have obtained that 


^-^[N{v)Y = ^ 
3 


i=l 

i=l 


-t d 1 -h a -h /3 D 

fll + a + 


a 


a 


i=l 


which exactly shows that ^ is paracontrolled by the collection 

[> 


Corollary 29. Under the assumptions of Theorem [2S1 on the positive parameters 
a,/3,a,e, and given u e Sa,j 3 {C) with u e set v := — {u ■ V)u). 

Then the tuple (y,u,ui,U 2 ) satisfies the structure equation (15.4p . with 


{v,U,Ul,U2) 


a,l3 




(u,Ui,U2,U3) 


a,/3 


(5.14) 


where k is the constant appearing in the definition (|5.10l) of the weight w. 

We summarise in the following assumption the work about renormalisation of the 
ill-defined terms defining the Burgers-enhancement of 

Assumption (B’) Assumption (B) hold and denoting by the quantity obtained by 
replacing ^ by in the definition (I5.12p of 0, then for any a < 1/2, there exists some 
deterministic 3x3 matrix-valued function such that 0^ — converges in probability 
in 


This assumption is the final ingredient in the proof of Theorem [2j 


Proof of Theorem [5] - Well-posedness of Burgers system follows as a direct con¬ 
sequence of Theorem [55J Theorem [5] on the convergence of the solutions to 
a renormalised e-dependent equation to the solution of the Burgers equation 
is thus obtained as a direct consequence of this well-posedness result together 
with an additional renormalisation step that will be done in a forthcoming work. 
The 3x3 matrix-valued functions d^ is the one renormalizing the quantities 
By tracking the changes (in the proof of Theorem [28P , induced by 
a renormalisation of 0 ^ into 0 ^ — # in we see that if (u^, uf, uli'^ 3 ) 

satisfies Ansatz (15.41) with Zf, and setting := • V)u^ — d^(uf,uf)), 

the tuple (u^, uf, ul) still satisfies the ansatz. We then complete the proof 
of Theorem [2l as done for Theorem [TJ 

[> 
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6 _ 

A glimpse at renormalisation matters 

We provide in this section a flavour of the problems that are involved in proving 
that assumption (B), formulated in Section 15 .31 hold true. The analysis of the 2 and 
3-linear terms is essentially complete, while the analysis of the 4-linear terms is only 
sketched. Hairer uncovered in [2U] the rich algebraic setting in which renormalisation 
takes place within his theory of regularity structures. The full treatment of this 
problem was given very recently in the works [10] and |15j of Hairer and co-authors. 
They provide in particular a clear understanding of which counterterms need to/can 
be added in the dynamics driven by a regularized noise to get a converging limit 
when the regularizing parameter tends to 0. A similar systematic treatment of 
renormalisation matters within the setting of high order paracontrolled calculus |5| 
should be developed in a near future. We describe in this last section how things 
can be understood from a pedestrian point of view on the example of the (PAM) 
equation. We assume here as in Theorem [1] that the vector fields V) are divergence- 
free-, this specific assumption is used to see that the terms do not need to be 
renormalised. 

Basic renormalisation consists in removing from diverging random terms their 
expectation. While this operation is sufficient in a number of cases, such as the 2 
and 3-dimensional (PAM) equations, or the 1-dimensional stochastic heat equation 
|26[ I28j . more elaborate renormalisation procedures are needed in other examples, 
such as the (KPZ) or $3 equations. Hopefully, the kind of renormalisation needed 
here for the study of the 3-dimensional (PAM) and Burgers equations, is essentially 
basic, in accordance with the work of Hairer and Labbe [28] on the (PAM) equation 
in 


The a priori ill-defined terms are 2-linear with respect to the noise 

and (vfiZi) , Zi)^, 

3-linear 

U^^'>{V^ZuViZ2), 

and 4-linear 

U^^\V,Zl,V^Zs), and H^'’) (H^Zs, P^Za). 


6.1. Renormalising the quadratic terms One takes advantage in the analysis of the 

renormalisation of the quadratic terms of 
the fact proved along the proof of Proposition [15] that the operator t~^ Qf is also 
a Gaussian operator with cancellation, an element of GCs” actually. More generally, 
the operators Qt o p o are of the form \ft Qfi for some Gaussian operator Q[ 
with cancellation. Thus the term H*-^^ (pZi, pZi) has the same structure as 

\ 2 --= 

so does the resonent term H^^^ (f, Zi) analysed in [3] in the study of the 2 -dimensional 
(PAM) equation. We estimate the size of Qrih) in terms of r, to see whether or not 
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it belongs to some parabolic Holder space. For a space white noise C, the expectation 
IE |Qr(l 2 )(e)|^ is given by the integral on x [0,1]^ of 

Xs,P,^(e,eOi^S.n,(e,e'OE[Q,\C(eOQU(eOQU(e'OQU(e")l (6-1) 

against the measure v{de')v{de")dtidt 2 - The expectation in (|6.ip is estimated with 
Wick’s formula by 

+ E[Q,\e(eOQ?,C(e")]E[QU(e")Q?iC(e')] 

^ (^ 1 ^ 2 ) + ^^l+^ 2 (e^ 

where d is the homogeneous dimension of the ambiant space M; the term 

comes from the first product of expectations. The quantity E |Qr(l 2 )|^ can thus 

be bounded above by the sum of two integrals, with dm := v[de')v[de")dtidt 2 and 
a first integral equal to 


J ^QrVt^ (e> e')KQ,Vt 2 (e, e")^q+42(e', e"fdm 


< 


< 


\ 


ti 


r + t 2 


(fl +^ 2 ) 2 J gr+ti{x,y)Gr+t2{x,z)Gti+t2{y,z)dm 


/(^) (^) it^ + hrHr + h + t.r^/^dhdt, 


^ r 


2-d 


for d < 4; we used the upper bound (13.41) here. We also have 

J J KQ^'p,^{e,e')KQ^'p,^{e,e"){tit2)~'^/‘^dv{de')v{de")dtidt2 


< 


I 




h 


{tit2) '2dtidt2 


for some positive exponent a, with a relatively sharp upper bound, which happens to 
be infinite in dimension 2 or larger. This is the annoying bit. Considering I 2 — E[l 2 ] 
instead of I 2 removes precisely this diverging part in the corresponding Wick formula 
|2l 


for E 


|Qr-(l 2 -E[l 2 ])| 


. It follows as a consequence that one has 


E 


|Qr-(l 2 -E[l 2 ])| 


<C 7" 


1 -# 


which shows that the associated distribution is almost surely in , by Kol¬ 

mogorov’s continuity criterion. 


While the above reasoning shows that recentering 

^0 

2=1 

around its expectation makes it converge in the right space, there is actually no need 
to renormalize this term, as can be expected from comparing our setting with the 
setting of regularity structures for the 3-dimensional setting, investigated in Hairer 
and Labbe’s work [28]. 
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One can see that point by proceeding as follows. Replace in a first step the study 
of the above quantity by a similar quantity where the spacetime paraproduct 
and resonent term are replaced by a space paraproduct and resonent 

operator •) introduced and studied in [1] - they are defined in the exact same 

way as II^, but without the time convolution operation. Continuity properties were 
proved for such spatial paraproduct in [3] , and we shall use in addition an elementary 
comparison result between this spatial paraproduct and our space-time paraproduct 
proved by Gubinelli, Imkeller and Perkowski in their setting |22l Lemma 5.1]. A 
similar statement and proof holds with the two paraproducts and ; we state 
it here for convenience. 

Lemma. Let uji,uj 2 be two space-time weights. If u e for a e (0,1) and v e L‘^Cu 2 
for some /3 e (—3,3) then 

7rW(u) - n(f)(u) e 


with w = u;iu; 2 . 

Setting 

i=l 

and using the comparison lemma and then the continuity estimates of each para¬ 
product, we see that is equal to 




2 = 1 






GO/^2a —1 


2=1 


so it is an element of So in order to estimate is the suitable Holder 

space we only need to study its ’’spatial” counterpart This can be done as 

follows. 


As the quantity is quadratic as a function of the noise, however we are 

going to see that its expectation is already bounded in as a consequence of 

some symmetry properties - this explains why is directly converging in 
with no renormalisation needed along the way. The term can indeed be written 
as a finite sum of integrals in time of terms of the form 

Pt (e) + Pt (e), 

where the localizing operators Pt and Qt are only in space. Using the above addi¬ 
tional geometric assumptions on the operator, the previous integral can be estimated, 
up to a satisfying remainder term controlled in terms of by 

Pt (e) + Pt (e). 

Its expectation can be seen to converge in to 

J Apt(a;,y) |^A[y,Q2^-i(y,.sf-i)]*Qjyi^-i(y,2/) + K[VjQl^-^{V3^-^)]*QlViS?-^iyTy)\ K^y), 
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where * denotes the usual adjoint in Lp‘[M,d^) (in space) and where the time is 
fixed in the operator . By symmetry, it is equal to 



(y, y) 


Kdy) 


where Qt '■= Q^’*VjQl + Ql’*VjQ^ is antisymmetric. Since at time fixed, the spatial 
operator is self-adjoint, we deduce that is antisym¬ 

metric in space and so its kernel is vanishing on the diagonal. This shows as a 
consequence that E[u; 2 ^] is bounded in the parabolic Holder space 


6.2. Higher order terms The analysis of the 3-linear term H^^^ (HjZi, 14 ^ 2 ) can be 

done exactly as for the 2 -linear term H^^^ {ViZi, VjZi), start¬ 
ing from the fact that the former has the same structure as 


\- 


■■= Vt, (qU • Ql {VtAQU • QU})) dt2dh, 


where the Ql are Gaussian operators with cancellation. Its renormalised version I 3 
is defined by replacing by ■ Q^C - ^[QtaC • Q^C]- The quantity 

r I 121 

E |Qr(l 3 )(e)| is thus given by an integral with respect to some kernels with some 
Gaussian controls and cancellation property of the expectation of a product of six 
Gaussian random variables indexed by parabolic points {e',e",e',e") 

{QlOie') (Q.\C)(e") {QUW) (QtC)(eO {QlC){e") {QUW')- 

In its renormalised version, the above product is replaced by 

{QU)(e') (0,C)(g') - E[(C?,C){e') 

and similarly for (Qs 2 C)(e ^0 Using Wick formula then shows that 

| 2 ‘ 


E 


Q.(IS)(e) 


only involves products of expectations where no two identical parameters appear 
inside each expectation, meaning that we have, after some elementary computations, 
an estimate of the form 

with dm = iy{de')i'{de")dtidt 2 , as above. For d < 4, this gives the estimate 

iE [|Qr(l3)(e)P] ~ J J + + h + t2)~'^^‘^tit2dtidt2 

< ^-3d/2+4^ 


on which one reads that I 3 has almost surely regularity (—3d/2 -I- 4)so no renor¬ 
malisation is required. This is coherent with [28] - equation (5.3), where it is shown, 
within the setting of regularity structures, that the terms that are trilinear functions 
of the noise do not need to be renormalised. Everything happens here as if we were 
working with a 3-linear term of the form 

frt{Qk-QU-tQlC)df, 

Jo 

one can indeed make that comparison concrete. 
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• The model quantities corresponding to the 4-linear terms are of the type 

14:= frt{Qk-Qk-tQk-tQtC)dt. 

Jo 

One can see on such terms that a basic renormalisation procedure suffices to get 
objects of regularity 0~, in dimension 3, such as expected. This finishes the sketch 
of proof that assumption (B) actually holds true. 
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